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Equations  for  a uniform  transmission  line  with  frequency-dependent  losses  were 
used  to  predict  the  performance  limits  for  an  optimized  transducer.  To  accountl 
for  transducer  non-uniformities,  we  modeled  the  line  as  a number  of  short 
segments  with  equal  lengths  but  different  impedances.  We  analyzed  the  loosely- 
laid  wire-in-tube  configuration  as  well  as  an  optimized  geometry  with  random 
variations  in  the  conductor  spacing. 

It  is  theoretically  possible  to  achieve  a range  resolution  of  30  meters  and  a 
signal-to-noise  ratio  of  25  dB  with  a 1/2  inch  diameter,  2 km  long  transducer. 
Range  resolution,  as  well  as  signal  amplitude,  are  limited  by  attenuation  and 
multiple  reflections.  Multiple  reflections  also  produce  clutter  and  cross 
talk,  which  limits  the  sensor  dynamic  range.  These  problems  become  quite  signi- 
ficant if  the  variations  in  conductor  spacing  exceed  about  one  millimeter  peak- 
to-peak  over  a spatial  scale  of  one  or  two  meters. 

Because  of  multiple-reflection  problems,  the  loosely-laid  wire-in-tube  trans- 
ducer is  limited  to  usable  lengths  of  a few  hundred  meters  in  pulsed-mode 
operation.  We  have  identified  a transducer  geometry  in  which  the  inner  con- 
ductor is  supported  within  an  eccentric  hollow  dielectric.  This  type  of  sensor 
should  alleviate  the  multiple  reflection  effects. 
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TECHNICAL  REPORT  SUMMARY 


PULSED-MODE,  ENCLOSED-WIRE  TRANSDUCER  ANALYSIS 


The  pulsed-mode,  enclosed-wire  transducer  offers  the  potential  for  signif- 
icant improvement  over  existing  buried-line  sensors.  Other  sensors  must 
be  deployed  in  multiple  segments  to  reduce  the  cumulative  effects  of  back- 
ground noise  and  to  provide  information  on  the  location  of  an  intrusion. 
Operating  in  the  pulsed  mode,  time  gating  can  be  used  to  separate  signals 
from  different  portions  of  the  transducer.  A single  long  transducer  and 
one  set  of  processing  electronics  can  provide  multiple  range-channel  out- 
puts. 

Analysis  and  optimization  of  this  sensor  concept  involves  a number  of 
parameters,  and  can  best  be  handled  by  computer  simulation.  This  study 
was  intended  to  provide  a basis  for  the  subsequent  design  and  feasibility 
testing  of  a breadboard  sensor. 

We  applied  transmission  line  theory  to  the  pulsed-mode  transducer  analysis. 
Although  the  work  was  primarily  analytical,  we  performed  limited  measure- 
ments on  an  experimental  transducer  section  to  help  in  establishing  a real- 
istic computer  simulation.  The  transducer  analysis  was  handled  in  two 
parts.  We  used  equations  for  a uniform  transmission  line  with  frequency- 
dependent  losses  to  predict  the  performance  limits  of  an  optimized  trans- 
ducer. For  the  case  where  the  conductor  spacing  is  not  uniform,  we  had 
to  develop  a computation  technique  w’hich  minimizes  the  bookkeeping 
involved  in  calculating  all  of  the  multiple  reflections.  The  technique  we 


used  treats  the  transducer  as  a lossless  transmission  line  made  up  of  many- 
short  segments  with  equal  lengths  but  different  impedances.  We  analyzed 
the  loosely-laid  wire-in-tube  configuration  as  well  as  an  optimized  trans- 
ducer geometry  with  random  variations  in  the  conductor  spacing. 

Our  results  can  be  briefly  summarized  as  follows: 

• It  is  theoretically  possible  to  achieve  range  resolutions  of 

30  meters  and  signal-to-noise  ratios  of  25  dB  with  a 1 /2-inch 
diameter,  2 km  long  transducer  and  a simple  processor. 

• Maximum  usable  transducer  length  is  limited  by  attenuation 
and  multiple  reflections.  These  factors  affect  the  range 
resolution  as  well  as  the  signal  amplitude.  Multiple  reflec- 
tions due  to  variations  in  conductor  spacing  can  also  cause 
crosstalk  between  range  channels. 

• The  outer  conductor  should  be  as  large  as  possible  to  mini- 
mize line  attenuation.  However,  such  factors  as  cost, 
weight  and  difficulty  of  deployment  place  a practical  size 
limit  on  the  order  of  one  inch  on  the  transducer  diameter. 

• Optimum  electrical  performance  is  obtained  with  an  inner 
conductor  diameter  which  is  about  0.  3 times  the  outer  con- 
ductor diameter.  Performance  is  relatively  constant  for 
diameter  ratios  between  0.  15  and  0.  6. 

• A gap  of  approximately  one  millimeter  between  the  inner 
and  outer  conductors  represents  an  acceptable  compromise 
between  sensitivity  and  resolution. 


• Variations  in  the  gap  must  be  minimized,  both  in  magnitude 
and  spatial  frequency.  When  the  gap  varia^ons  exceed  about 
one  millimeter  peak-to-peak  over  a spatial  scale  of  less  than 
one  or  two  meters,  clutter  and  crosstalk  become  significant 
problems. 

• Because  of  clutter,  crosstalk  and  attenuation  caused  by 
multiple  reflections,  the  loosely-laid  wire-in-tube  transducer 
is  not  well  suited  for  pulsed-mode  operation.  This  trans- 
ducer geometry  is  unlikely  to  be  usable  in  lengths  greater 
than  a few  hundred  meters.  How-ever  there  are  alternate 
transducer  configurations  w'hich  appear  to  satisfy  the  con- 
flicting design  requirements --simplicity  and  ruggedness,  low- 
cost,  uniform  conductor  spacing,  and  sensitivity  to  mechani- 
cal disturbances. 

The  next  logical  step  in  the  development  of  the  pulsed-mode,  enclosed-wire 
transducer  is  an  experimental  confirmation  of  the  analytical  predictions. 

One  or  more  candidate  transducer  configurations  should  be  fabricated  and 
then  tested  under  actual  field  conditions.  These  tests  should  confirm  the 
feasibility  of  the  pulsed-mode  transducer  concept.  A combination  of  experi- 
mental and  analytical  data  will  provide  a firm  basis  for  further  refinements 
in  the  transducer  and  processor  design. 
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EVALUATION 


This  report  provides  valuable  insight  into  the  maximum  feasible  length 
transducer  to  which  the  pulsed  mode  concept  may  be  applied.  The  contractor 
identified  early  in  the  study  the  inherent  limitations  associated  with  pulsing 
a wire  loosely-laid  inside  of  a copper  tube.  Severe  cross  talk  dispersion 
and  attenuation  plagued  this  concept,  and  a minimum  length  of  only  a few 
hundred  meters  could  be  anticipated.  Honeywell  then  showed  how  to  optimize 
the  geometry  of  the  enclosed  wire  transducer  and  calculated  the  performance 
limitation  for  transducers  with  diameters  of  one  inch.  Results  of  the  study 
of  an  offset  derandoraized  version  of  the  line  indicate  that  a 100  meter 
penetration  location  capability  can  be  achieved  from  transducer  which  can 
extend  some  3 km  in  length.  The  information  provided  by  the  report  will  be 

I 

used  in  subsequent  design  and  feasibility  testing  of  a breadboard  sensor. 

This  effort  has  been  funded  entirely  through  Laboratory  Directors'  Funds 

uiuler  TP0-R6A  in  support  of  the  Base  and  Installation  Security  System  Program. 
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SAMUEL  J.  STECHER,  Jr.  Lieutenant,  USAF 
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viii 


SECTION  I 


INTRODUCTION 


Honeywell  has  completed  an  analysis  of  the  Pulsed-Mode,  Enclosed-Wire 
Transducer  for  Rome  Air  Development  Center.  This  transducer  concept 
offers  the  potential  for  significant  improvement  over  existing  buried-line 
sensors.  Other  sensors  must  be  deployed  in  multiple  segments  to  limit 
the  effects  of  background  noise  and  to  provide  information  on  the  location  of 
an  intrusion.  Operating  the  enclosed-wire  transducer  in  the  pulsed  mode, 
time  gating  can  be  used  to  separate  signals  from  different  portions  of  the 
transducer.  A single  long  transducer  and  one  set  of  processing  electronics 
can  provide  multiple  range-channel  outputs. 

The  objective  of  our  study  was  to  provide  a basis  for  the  subsequent  design 
and  feasibility  testing  of  a breadboard  sensor.  This  report  shows  how  to 
optimize  the  transducer  design  for  pulsed  operation  and  predicts  the  per- 
formance attainable  with  practical  transducer  configurations.  All  of  the 
work  was  performed  at  Honeywell's  Systems  and  Research  Center. 
Minneapolis,  under  the  direction  of  Dr.  George  Webber. 

In  our  study,  we  applied  transmission-line  theory  to  the  analysis  of  the 
enclosed-wire  transducer  operating  in  the  pulsed  mode.  The  work  was 
primarily  analytical,  although  we  measured  the  characteristics  of  an 
experimental  transducer  to  form  a basis  for  the  computer  simulation.  The 
transmission-line  analysis  was  handled  in  two  parts.  We  analyzed  a uni- 
form transmission  line  with  frequency-dependent  losses  to  predict  the 
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performance  limits  of  an  optimized  transducer.  For  the  case  where  the 
conductor  spacing  is  not  uniform,  we  had  to  develop  a computation  technique 
which  minimizes  the  bookkeeping  involved  in  calculating  all  of  the  multiple 
reflections.  The  technique  that  we  used  treats  the  transducer  as  a lossless 
transmission  line  made  up  of  many  short  segments  with  equal  lengths  but 
different  impedances.  We  analyzed  the  loosely-laid  wire-in-tube  configur- 
ation as  well  as  an  optimized  transducer  geometry  with  random  variations 
in  the  conductor  spacing. 

Our  results  show  that  an  optimized  ^-inch  diameter,  2-km  long  transducer 
with  a simple  processor  can  meet  or  exceed  Air  Force  performance  require- 
ments. However,  the  loosely-laid  wire-in-tube  configuration  suffers  a 
severe  performance  degradation  because  of  multiple  reflections.  We  have 
identified  alternate  transducer  designs  which  should  alleviate  this  problem. 


Included  in  this  report  are  results  of  our  transducer  analysis  and  design 
optimization  studies.  There  is  a brief  discussion  of  signal  processing  with 
a sample  signal-to-noise  calculation.  Several  alternate  transducer  designs 
are  recommended  for  a possible  breadboard  feasibility  verification. 


SECTION  II 


SUMMARY  AND  CONCLUSIONS 


During  this  study,  we  showed  how  to  optimize  the  geometry  of  the  pulsed- 
mode,  enclosed-wire  transducer,  and  we  calculated  the  performance  limits 
for  transducers  with  diameters  up  to  one  inch.  We  also  analyzed  the  effects 
of  non-uniformities  in  the  conductor  spacing.  Our  results  can  be  briefly 
summarized  as  follows: 

• It  is  theoretically  possible  to  achieve  range  resolutions  of  30 
meters  and  signal-to-noise  ratios  of  25  dB  with  a 1/2  inch  diame- 
ter, 2 km  long  transducer  and  a simple  processor. 

• Maximum  usable  transducer  length  is  limited  by  attenuation  and 
multiple  reflections.  These  factors  affect  the  range  resolution  as 
well  as  the  signal  amplitude.  Multiple  reflections  due  to  vari- 
ations in  conductor  spacing  can  also  cause  crosstalk  between 
range  channels. 

• The  outer  conductor  should  be  as  large  as  possible  to  minimize 
line  attenuation.  However,  such  factors  as  cost,  weight,  and 
difficulty  of  deployment  place  a practical  size  limit  on  the  order 
of  one  inch  on  the  transducer  diameter. 

• Optimum  electrical  performance  is  obtained  with  an  inner  conduc- 
tor diameter  which  is  about  0.  3 times  the  outer  conductor  diame- 
ter. Performance  is  relatively  constant  for  diameter  ratios  be- 
tween 0.15  and  0.6. 
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• A gap  of  approximately  one  millimeter  between  the  inner  and  outer 
conductors  represents  an  acceptable  compromise  between  sensi- 
tivity and  resolution. 

• Variations  in  the  gap  must  be  minimized,  both  in  magnitude  and 
spatial  frequency.  When  the  gap  variations  exceed  about  one  milli- 
meter peak-topeak  over  a spatial  scale  of  less  than  one  or  two 
meters,  clutter  and  crosstalk  become  significant  problems. 

• Because  of  clutter,  crosstalk  and  attenuation  caused  by  multiple 
reflections,  the  loosely-laid  wire- in- tube  transducer  is  not  well 
suited  for  pulsed-mode  operation.  This  transducer  geometry  is 
unlikely  to  be  usable  in  lengths  greater  than  a few  hundred  meters. 
However,  there  are  alternate  transducer  configurations  which 
appear  to  satisfy  the  conflicting  design  requirements:  simplicity 
and  ruggedness,  low  cost,  uniform  conductor  spacing,  and  sensi- 
tivity to  mechanical  disturbances. 

The  next  logical  step  in  the  development  of  the  pulsed-mode,  enclosed-wire 
transducer  is  an  experimental  confirmation  of  the  analytical  predictions. 

One  or  more  candidate  transducer  configurations  should  be  fabricated  and 
then  tested  under  actual  field  conditions.  These  tests  should  confirm  the 
feasibility  of  the  pulsed-mode  transducer  concept.  A combination  of  experi- 
mental and  analytical  data  will  provide  a firm  basis  for  refinements  in  the 
transducer  and  processor  design. 


SECTION  III 


TECHNICAL  DISCUSSION 


SENSOR  DESCRIPTION 


The  pulsed-mode,  enclosed-wire  transducer  uses  time-domain  reflecto- 
meter  techniques  to  provide  multiple  range  channel  outputs  from  a single 
long  sensor.  In  its  simplest  form,  the  transducer  can  be  an  insulated  wire 
loosely  laid  within  a conducting  tube.  Reference  1 describes  the  operation 
of  such  a sensor  under  d-c  bias  conditions.  The  capacitance  between  the 
wire  and  tube  is  a function  of  their  relative  positions.  Because  the  wire  is 
free  to  move  within  the  outer  tube,  a mechanical  disturbance  can  cause  a 
relative  motion  of  the  conductors  and  a corresponding  time  variation  in  the 
transducer  capacitance.  If  the  transducer  is  biased  with  a d-c  voltage 
through  a large  resistor,  the  capacitance  change  results  in  a time-varying 
voltage  which  can  be  detected  with  sensitive  electronics.  In  fact,  there 
may  be  enough  electret  charge  trapped  in  the  insulation  so  that  no  external 
biasing  is  necessary.  The  sensor  is,  in  effect,  a distributed  capacitance 
microphone  which  can  sense  disturbances  along  its  entire  length. 


Because  of  its  simplicity,  sensitivity,  and  low  cost,  the  enclosed-wire 
transducer  is  attractive  as  a buried-line  sensor  for  perimeter  defense. 
However,  several  factors  limit  the  maximum  transducer  length  to  some- 
thing on  the  order  of  100  meters.  Sensitivity  is  reduced  as  the  line  length 
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is  increased,  because  of  the  increase  in  total  shunt  capacitance.  Cumula- 
tive effects  of  background  noise  in  a long  distributed  sensor  eventually 
swamp  out  the  signal.  Furthermore,  it  is  usually  desirable  to  know  the 
location  of  an  intrusion  to  within  100  meters  or  less. 

If  pulsed  excitation  is  used  rather  than  a d-c  bias,  intrusion  signals  from 
various  points  along  the  transducer  will  be  delayed  by  different  amounts  of 
time  following  the  excitation  pulse.  In  this  mode  of  operation,  the  sensor 
is  essentially  a closed-loop  radar  system.  Time  gating  can  be  used  to 
separate  the  signals  into  range  channels  for  additional  processing.  This 
makes  it  possible  to  suppress  background  noise  and  to  provide  many  range- 
channel  outputs  from  one  continuous  long  transducer  and  a single  set  of 
processing  electronics. 

THE  ENCLOSED-WIRE  TRANSDUCER 
AS  A TRANSMISSION  LINE 

Basically,  the  enclosed-wire  transducer  operating  in  the  pulsed  mode  can 
be  treated  as  a non-ideal  transmission  line.  Using  equations  which  relate 
the  transducer  geometry  to  the  characteristic  line  impedance,  transmission- 
line  theory  can  be  applied  directly  to  a prediction  of  transducer  perfor- 
mance. 

II 

Wave  Propagation  on  a Transmission  Line 

if  a signal  is  applied  to  a transmission  line,  it  propagates  down  the  line  at 
a finite  velocity  so  that  the  voltage  at  point  x lags  the  phase  of  the  input 
voltage  by  an  amount  9 per  unit  length.  The  amplitude  of  the  wave  at  point 
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x in  a non-ideal  line  is  also  reduced  by  an  amount  a per  unit  length  because 
of  conductor  resistance  and  dielectric  losses.  A complex  propagation  con- 
stant v is  defined  as: 


V = or  + j 9 


where: 

j = -£l 


At  any  point  on  a uniform,  infinitely  long  line  the  voltage  and  current  are 
given  by: 


E(x)  = E.  e‘VX  I(x)  = I.  e'VX 
in  in 

A transmission  line  has  a characteristic  impedance  Z which  is  determined 

o 

by  the  line  geometry  and  the  properties  of  the  conductors  and  dielectric 
material  in  the  line.  On  an  infinitely  long  line,  the  ratio  of  voltage  to  cur- 
rent is  the  same  at  all  points  and  is  equal  to  Z^.  This  is  also  true  if  the 
line  is  finite  in  length  and  is  terminated  in  a lumped  impedance  equal  to  Zq. 
But  if  the  line  is  not  "matched",  some  energy  is  not  delivered  to  the  load 
and  is  accounted  for  by  a reflected  wave  that  propagates  back  toward  the 
source.  This  reflection  condition  will  also  exist  if  the  line's  characteristic 
impedance  is  not  constant  along  its  path.  Figure  1 illustrates  a voltage  step 
propagating  from  a source  down  a transmission  line,  and  the  corresponding 
reflected  voltage  from  a discontinuity  in  the  line  impedance. 


CHARACTERISTIC  IMPEDANCE  CHARACTERISTIC  IMPEDANCE 


Figure  1.  Reflection  from  a Discontinuity  in  Line  Impedance 

If  the  input  pulse  is  applied  to  the  line  at  time  tQ,  the  reflected  wave  from  an 
impedance  discontinuity  at  point  x arrives  back  at  the  beginning  of  the  line  at 
a time  t = t + 2x/c:  where  c is  the  propagation  velocity  of  an  electromag- 
netic wave  on  the  line.  In  an  air-insulated  transmission  line,  c is  very  close 

g 

to  the  speed  of  light  in  a vacuum,  3 x 10  meters  per  second. 
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The  ratio  of  the  reflected  voltage  to  the  incident  voltage  V.,  called  the 
voltage  reflection  coefficient  p,  is  related  to  the  change  inline  impedance  by 

Z ' - Z V 
_ _o o r 

D “ Z ' + Z ~ V. 
o o 1 

N'ote  that  if  Z 'is  less  than  Z , p is  negative,  indicating  that  the  wave  re- 
o o 

verses  sign  upon  reflection.  This  will  be  important  later  in  considering  the 
reflection  from  a small  disturbance  in  the  line  impedance.  We  also  note 
that  for  small  changes  in  line  impedance,  the  reflection  coefficient  can  be 
approximated  by 

~ AZ 

A ■ 

” 2Z 

o 

Impedance  of  a Transducer  with  a Uniform  Dielectric 

As  mentioned  above,  the  characteristic  impedance  of  a transmission  line  is 
determined  by  the  goemetrv  of  the  line  and  the  properties  of  the  conductors 
and  dielectric.  Some  form  of  dielectric  coating  is  needed  on  the  inner  con- 
ductor to  prevent  a short  in  the  line  if  a sharp  bend  is  introduced.  If  the 
dielectric  layer  is  fairly  thin,  it  will  have  little  effect  on  line  impedance  or 
sensitivity.  A reasonably  good  approximation  to  the  line  impedance  can  be 
made  by  assuming  a uniform  dielectric  constant.  Computation  procedures 
for  the  impedance  of  an  eccentric  line  with  an  inhomogeneous  dielectric  are 
discussed  in  Appendix  B. 

The  resistivity  of  the  conductors  is  typically  low  enough  to  be  neglected  in 
the  calculation  of  characteristic  impedance,  even  though  it  gives  rise  to  an 
attenuation  factor  which  can  by  no  means  be  neglected.  For  the  eccentric 
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transducer  geometry  the  characteristic  impedance  is  not  only  a function  of 
the  ratio  of  the  outer  conductor  to  inner  conductor  diameter,  but  also  depends 
on  the  distance  of  the  inner  conductor  from  the  axis  of  the  outer  tube  (Figure 
2). 


Figure  2.  Transducer  Geometry 

In  the  determination  of  characteristic  impedance,  we  will  assume  that  the 
conductors  and  dielectric  have  the  same  permeability  as  free  space.  The 
impedance  of  a transmission  line  can  be  expressed  in  terms  of  its  inductance 
per  unit  length  L and  its  capacitance  per  unit  length  C by 
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The  velocity  of  propagation  in  the  transverse  electromagnetic  mode  is 


V = i 

(e  )? 
r 

where  c is  the  velocity  of  tight  in  free  space,  and  is  the  relative  permit- 
tivity of  the  dielectric  material  filling  the  transmission  line.  The  velocity 
can  also  be  expressed  as 


V = (LC) 


so  that  we  can  define  Z as 

o 


Z = 


(e  )■ 
r 


° cC 


The  capacitance  per  unit  length  of  an  eccentric  cable  uniformly  filled  with  a 
dielectric  material  is 


C = 


2 TT  e e 
o r 


cosh  1 [(aj “ + a9"  - d')/2a  a,)] 


-1° 

where  eQ  = 8.  85  x 10  “ farads  per  meter.  Combining  the  expressions,  we 

can  write  Z as 
o 


Z = 60(e  ) ‘ cosh  *[(a“  +a  “ - d“)/2a,a„] 
o r 12  12 
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versus  inner-conductor  position  for  an  air-dielectric  cable  consisting  of  a 
0.  14-inch  diameter  wire  in  a f-inch  ID  tube.  Also  shown  is  the  reflection 
coefficient  resulting  from  a 0.01  millimeter  change  in  the  inner-conductor 

1 

position.  We  have  expressed  the  position  in  terms  of  the  gap  g between  the 

conductors,  which  is  related  to  the  other  dimensions  by  . j 1 

d = a2  " (al  + g) 


The  sensitivity  to  small  relative  displacements  is  zero  when  the  inner  con- 
ductor is  centered,  and  reaches  a maximum  when  the  inner  conductor  is 
nearly  touching  the  tube  wall. 

Signal  Attenuation  in  the  Enclosed-Wire  Transducer 


Even  in  an  air-insulated  line  with  no  dielectric  losseo,  the  finite  conductivity 
of  the  conductors  results  in  a loss  of  amplitude  as  a wave  prppagates  down 
the  line.  Furthermore,  the  losses  become  greater  at  higher  frequencies. 
This  is  due  primarily  to  the  "skin  effect"  which  confines  the  high-frequency 
currents  to  the  surface  region  of  a conductor  and  gives  rise  to  a frequency- 
dependent  a-c  resistance.  Although  the  exciting  waveform  is  a d-c  voltage 
change,  we  cannot  assume  that  only  d-c  losses  are  important.  Actually,  a 
Fourier  analysis  of  an  instantaneous  step  voltage  change  reveals  that  it  con- 
tains an  infinite  band  of  frequencies.  If  the  step  is  applied  to  a transmission 
line  with  frequency-dependent  losses,  the  rise  time  of  the  step  waveform 
will  continuously  increase  as  it  propagates  down  the  line.  As  we  will  show 
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Figure  3.  Characteristic  Impedance  and  Reflection  Coefficient  versus  Gap 
in  a 1 /2-Inch  Transducer 


later,  a degradation  in  rise  time  results  in  a reduction  in  the  amplitude  of 
the  reflection  from  a small  disturbance  in  the  line.  It  also  "smears"  the 
reflected  signal  out  over  a long  time  period,  limiting  the  spatial  resolution 
of  disturbances  on  the  line. 

.A  calculation  of  the  attenuation  at  a single  frequency  does  not  yield  direct 
information  on  the  signal  amplitudes  and  waveforms  in  the  pulsed-mode 
transducer.  However,  the  equations  describing  the  attenuation  provide  in- 
sight into  the  parameters  which  affect  the  transducer  response. 

Factors  which  influence  the  attenuation  in  the  line  include  the  conductivity  of 
the  inner  wire  and  the  shield,  and  the  diameter  of  the  conductors.  Since 
current  density  is  higher  in  the  inner  wire  than  in  the  shield,  its  resistive 
losses  usually  dominate  in  determining  the  line  attenuation.  Losses  also 
vary  inversely  with  the  characteristic  impedance  of  the  line.  The  equations 
which  follow  are  for  an  air-insulated  line.  They  are  also  reasonably  accur- 
ate if  the  line  is  filled  with  a low-loss  dielectric. 

Attenuation  per  unit  length  a in  a transmission  line  is  given  by 


where  R is  the  resistance  per  unit  length  at  the  frequency  of  interest.  This 
is  related  to  the  skin-effect  resistance  R by 


and  R^  at  a frequency  f is 


JLt  , 


F = (— )^ 

S CT 
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For  copper  conductors,  the  permeability  u = u = 4 TT  x 10  henries /meter 

7 ° 

and  the  conductivity  o in  5.  8 x 10  mhos/meter.  The  attenuation  for  a cop- 
per transmission  line  is  then 


2.  08  x 10  ,11  , 

y = = (f)  ( — + — ) nepers/meter 

O 12 


-‘■u  ^ ^ — (iV  (—  + — ) dB/kilometer 
o 12 


Figure  4 shows  the  attenuation  in  dB  per  km  versus  gap  for  a -|-inch  dia- 
meter transducer  at  a frequency  of  5 MHz.  This  is  approximately  the  band- 
width needed  for  a range  resolution  of  30  meters.  Since  the  drive  signal 
and  reflected  signal  are  both  attenuated,  the  total  effective  attenuation  per 
kilometer  of  transducer  length  is  twice  the  value  show-n.  Comparing  this 
curve  with  Figure  3,  we  note  that  the  high  sensitivity  at  small  gaps  is 
accompanied  by  high  attenuation.  For  a resolution  of  100  meters,  the  band- 
width can  be  reduced  by  a factor  of  3.  33  and  the  attenuation  would  be  re- 
duced by  the  factor  (3.  33)£  = 1.  82. 


MECHANICAL  RESPONSE  OF  THE 
ENCLOSED- WIRE  TRANSDUCER 


To  perform  a complete  analysis  of  the  pulsed-mode,  enclosed-wire  trans- 
ducer, it  is  necessary  to  know  the  exact  transmission-line  geometry.  We 
are  interested  in  the  static  variations  in  the  dimensions  of  a deployed  trans- 
ducer, as  well  as  the  conductor  motion  caused  by  a time-varying  disturbance. 
In  subsequent  discussions,  we  use  experimental  data  as  well  as  computer 
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5 MHz 

1.78  MM  = 0.07  INCHES 
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simulations  to  establish  mechanical  models  for  the  static  dimension  vari- 
ations. We  use  these  models  to  show  the  effects  of  non-uniformity  on  trans- 
ducer performance. 

It  is  very  difficult  to  make  a quantitative  estimate  of  the  transducer's  time- 
varying  response.  .As  pointed  out  in  Reference  1,  there  are  several  modes 
of  response,  depending  upon  the  disturbance  frequency.  .At  very  low  fre- 
quencies, the  dominant  cause  of  relative  motion  between  the  wire  and  the 
tube  is  a deformation  of  the  tube  between  inner-conductor  support  points. 
Mechanical  resonances  in  the  inner  conductor  and  its  support  structure  give 
rise  to  the  so-called  "strumming"  mode.  For  the  looselv-laid  wire-in-tube 
sensor,  these  resonances  lie  in  the  20  to  40  Hz  range.  When  the  excitation 
frequency  is  high  enough,  the  inner  conductor  remains  relatively  stationary 
as  the  tube  vibrates.  In  all  of  these  modes,  the  relative  motion  is  stronglv 
influenced  by  the  inner  conductor  size  and  stiffness,  the  way  in  which  it  is 
supported,  and  the  surrounding  soil  conditions. 

Except  for  resonant  effects,  the  relative  motion  of  the  inner  and  outer  con- 
ductors cannot  exceed  the  deformation  of  the  soil  in  which  the  transducer  is 
buried.  Thus,  we  can  estimate  the  maximum  amplitude  of  the  time-varying 
conductor  displacement  by  calculating  the  soil  motion  under  loading  condi- 
tions representing  a walking  intruder. 

.Appendix  .A  gives  the  results  of  a quasi-static  solution  for  the  soil  displace- 
ment produced  by  a 200-pound  load  distributed  uniformly  over  an  8-inch 
radius  disk.  At  a depth  of  eight  inches  directly  below  the  load,  the  vertical 
soil  displacement  is  722  microinches  (approximately  0.02  millimeter);  the 
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maximum  displacement  decreases  to  451  microinches  (0.011  millimeter)  at 
a burial  depth  of  16  inches.  At  40  inches  (1  meter)  from  the  center  of  the 
load,  the  vertical  soil  displacement  is  nearly  constant  at  110  microinches 
(0.  0028  millimeter)  down  to  burial  depths  of  24  inches. 

We  will  show  later  that  the  static  variations  in  inner-conductor  relative 
displacement  must  be  minimized  to  reduce  attenuation  and  crosstalk.  It  is 
impractical  to  keep  the  outer  conductor  of  a buried  transducer  straight. 
Therefore  it  will  probably  be  necessary  to  support  the  inner  conductor  at 
intervals  not  much  greater  than  one  meter  in  order  to  restrict  its  maximum 
relative  displacement.  Based  upon  these  considerations  and  the  soil- 
displacement  computations,  we  have  assumed  that  the  time-varying  inner- 
conductor  displacement  has  a peak  amplitude  of  400  microinches  (0.01  milli 
meters)  and  an  effective  length  of  one  meter. 


EXPERIMENTAL  CHARACTERIZATION  OF 
THE  ENCLOSED-WIRE  TRANSDUCER 

Our  study  was  primarily  analytical,  but  we  performed  limited  field  experi- 
ments to  help  define  a simulation  analytical  model  for  the  transducer.  A 
10-meter  section  of  transducer  was  fabricated  and  installed  in  the  line  sen- 
sor test  facility  at  the  Honeywell  Ordnance  Proving  Ground.  The  transducer 
was  the  basic  wire-in-tube  sensor,  a number  26  Teflon-insulated  wire  in  a 
^ -inch  refrigeration  tube.  Time-domain  reflectometer  measurements  were 
made  to  determine  the  static  impedance  variations  along  the  length  of  the 
transducer.  Results  of  these  measurements  and  their  incorporation  into 
the  transducer  model  are  treated  in  detail  in  the  discussion  of  non-uniform 
transducer  effects. 
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We  also  measured  the  time-varying  reflection  coefficients  in  response  to  a 
known  soil  disturbance.  A cyclic  signal  generator  was  used  to  impart  a 
10-Hz,  100-pound  peak-to-peak  load  to  the  soil  surface.  The  generator  was 
placed  above  the  midpoint  of  the  transducer,  which  was  buried  about  eight 
inches  deep  in  sand. 

The  time-domain  reflectometer  controls  were  set  to  display  the  reflection 
coefficient  at  a single  point  on  the  transducer,  selectable  by  varying  the 
delay-time  potentiometer.  The  vertical  deflection  output  of  the  TDR  was 
recorded  on  a strip  chart  recorder,  along  with  the  outputs  of  geophones 
adjacent  to  the  transducer  and  directly  beneath  the  10-Hz  cyclic  signal  gen- 
erator. Recordings  of  the  TDR  output  were  made  at  time-delay  increments 
of  5 nanoseconds,  corresponding  to  position  increments  of  about  0.  6 meters 
along  the  line  transducer.  With  a peak-to-peak  soil  displacement  of  300 
microinches,  as  indicated  by  the  vertical-axis  geophone,  the  peak-to-peak 
variation  in  reflection  coefficient  near  the  signal  generator  was  about 

_3 

8.  6 x 10  . The  noise  level  with  this  instrumentation  was  fairly  high,  and 

it  was  difficult  to  distinguish  the  signal  from  the  noise  at  points  more  than 
two  meters  from  the  signal  generator.  In  fact,  except  for  one  measurement 
point,  the  signal  appeared  to  consist  mostly  of  the  second  harmonic  of  the 
excitation.  This  implies  a "strumming"  mode  excitation  of  the  wire  at  a 
harmonic  of  the  mechanical  input  frequency. 

We  also  measured  the  velocity  of  propagation  in  the  transducer  by  using  a 
short-circuit  load  at  the  far  end  to  provide  distinct  time  marker.  The 
round-trip  propagation  time  within  the  10-meter  test  section  was  80  nano- 
seconds. This  corresponds  to  an  average  propagation  velocity  of  0.  833 
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times  the  free-space  velocity  of  light.  The  average  velocity  factor  is  about 
midway  between  the  free-space  value  and  the  velocity  factor  of  a cable  with 
solid  Teflon  dielectric. 


With  the  transducer  shorted,  the  fall  time  of  the  reflected  step- function 
signal  provides  a measurement  of  the  round-trip  waveform  distortion  in  the 
transducer.  The  time  for  the  reflected  step  to  decay  to  50  percent  of  its 
final  value  was  approximately  three  nanoseconds.  This  is  much  longer  than 
the  fall  time  calculated  on  the  basis  of  frequency-dependent  line  losses. 


For  a step  excitation,  the  50  percent  risetime  ( 50%)  of  a coaxial  line  is 
approximately  given  by  the  formula 


T u 

50%  * 16  n-r2  ctZ  2 

o 


where  jj.  = permeability  of  the  inner  conductor 

l = length  of  the  cable 

r = radius  of  the  inner  conductor 

a - conductivity  of  the  inner  conductor 

Z = characteristic  impedance  of  the  cable 
o 


The  fall  time  of  the  reflected  step  on  a shorted  10-meter  line  is  obtained  by 
computing  the  rise  time  of  a step  excitation  on  a 20-meter  line.  For  a #26 


wire. 


p.  = 4TTx  10"7  H/M 
L = 20  m 

r = 0.  009  x 2.  54  x 10  m 
cr  = 5.  85  x 107  (0  meter)  * 
Z = 850 


T -q 

50%  = 0.  144  x 10  secs 

i 

This  calculation  led  us  to  conclude  that  the  pulse  degradation  on  the  test 
transducer  would  have  to  be  accounted  for  by  losses  due  to  multiple  reflec- 
tions. A valid  numerical  solution  of  the  line  response  must  therefore  show 
a fall  time  for  the  shorted-line  case  which  is  at  least  in  rough  agreement 
with  the  observed  value. 

TRANSDUCER  DESIGN  OPTIMIZATION 

From  the  standpoint  of  electrical  performance  of  the  sensor,  the  important 
quantities  are  the  amplitude  and  time  width  of  the  signal  resulting  from  a 
small  relative  disturbance  of  the  conductor  over  a short  section  of  its  length. 

A walking  intruder  may  produce  a typical  soil  displacement  of  about  0.01 
millimeters  at  a burial  depth  of  0.  25  meters  (10  inches),  and  the  disturbance 
extends  over  an  effective  distance  of  about  one  meter.  These  values  for  the 
transducer  displacement  were  used  in  our  design  optimization  studies. 
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Time  duration  of  the  signal  can  be  related  to  the  range  resolution  of  the 
sensor  by 


AR  = 

2 T 


where  c is  the  velocity  of  propagation  for  an  electrical  signal  in  the  trans 
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ducer  (for  air-insulated  line,  c = 3 x 10°  meters  per  second)  and  t is  the 
time  width  of  the  signal  from  a disturbance.  We  are  assuming  that  two 
pulses  can  be  resolved  if  their  separation  is  greater  than  their  width.  If 
the  transducer  had  no  losses  and  had  no  multiple  reflections,  the  width  of  a 
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signal  reflected  by  a disturbance  would  be  determined  only  by  the  length  of 
the  disturbance.  However,  the  effects  of  frequency-dependent  losses  and 
multiple  reflections  in  a real  transmission-line  transducer  not  only  decrease 
the  signal  amplitude,  but  also  broaden  it  in  time.  Transducer  geometry 
factors  influencing  the  signal  amplitude  and  duration  are  the  inner  and  outer 
conductor  dimensions,  the  nominal  gap  between  the  conductors,  the  amount 
of  variation  in  the  gap,  and  the  distance  from  the  generator  end  of  the  line 
to  the  disturbance. 


Signal  amplitude  in  the  pulsed-mode  transducer  is  directly  related  to  the 
drive  voltage.  However,  the  signal-to-noise  ratio  cannot  be  improved  in- 
definitely by  raising  the  drive  level.  A limit  is  reached  where  fluctuations 
in  the  drive  signal  become  the  dominant  noise  source.  Rather  than  express- 
ing signal  amplitude  in  absolute  units,  the  effective  reflection  coefficient  is 
a more  meaningful  term.  This  is  related  to  the  ratio  of  the  reflected  signal 
to  the  drive  signal. 

The  drive  voltage  level  used  in  our  computer  simulations  was  chosen  such 
that  the  signal  amplitude  is  numerically  equal  to  the  effective  reflection 
coefficient. 

Predicted  Signal  Waveforms 

Earlier,  we  showed  how  to  derive  the  characteristic  impedance  of  the  en- 
closed wire  transducer  as  a function  of  inner  conductor  position,  and  how 
the  voltage  reflection  coefficient  can  be  calculated  when  there  is  a step 
discontinuity  in  line  impedance.  Inner-conductor  displacement  with  respect 
to  the  outer  shield  will,  of  course,  be  gradual  rather  than  abrupt.  However, 


the  total  reflection  can  still  be  computed  with  the  same  equations  simply  by 
representing  the  inner-conductor  position  as  a series  of  small  discontinui- 
ties and  superposing  all  of  the  reflections  (with  the  appropriate  time  dis- 
placement) to  give  the  resultant  reflected  waveform.  Because  of  the  time 
"smearing"  of  the  incident  pulse  due  to  high-frequency  losses  in  the  line, 
information  on  the  exact  shape  of  a small  disturbance  in  the  center-conductor 
position  will  be  lost  on  a long  line.  In  fact,  it  is  a reasonably  valid  assump- 
tion to  approximate  the  gradual  center-conductor  position  change  by  a double 
step  as  shown  in  Figure  5.  This  reduces  the  computation  to  a superposition 
of  the  reflections  from  the  two  discontinuities. 


OUTER  SHIELD  ACTUAL  RELATIVE  POSITION 


DOUBLE-STEP  APPROXIMATION 
(DISPLACEMENT  = 0.01  mm) 


Figure  5.  Double-Step  Approximation  to 
Conductor  Displacement 
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Figure  5 shows  the  center  conductor  distorted  with  respect  to  the  outer 
conductor,  even  though  the  reverse  may  be  true.  Relative  displacement  is 
the  parameter  of  interest. 

Since  we  are  interested  in  the  reflected  wave  as  it  appears  back  at  the  input 
end  of  the  line,  the  attenuation  and  distortion  of  the  wave  in  propagating 
back  through  the  line  must  also  be  considered.  The  line  has  the  same 
propagation  constants  for  waves  traveling  in  either  direction,  so  we  can 
treat  the  reflection  from  the  discontinuity  at  point  L on  the  line  as  if  the 
incident  wave  at  the  location  of  the  discontinuity  is  simply  scaled  down  by 
the  linear  factor  p,  and  then  continues  along  an  identical  line  section  of 
length  L.  Mathematically,  it  makes  no  difference  whether  we  multiply  by 
p when  the  wave  gets  to  position  x or  wait  until  it  has  reached  2L.  Thus, 
the  reflected  waveform  appearing  at  the  driven  end  of  a pulsed  line  resulting 
from  a step  discontinuity  at  point  L is  equivalent  to  the  pulse  waveform  at 
point  2L  on  a uniform  line,  multiplied  by  the  voltage  reflection  coefficient 

of  the  discontinuity. 

With  the  double-step  approximation,  there  are  two  reflected  signals  separ- 
ated in  time  by  At  = 2A2/V.  For  small  displacements,  the  reflection  coeffic- 
ients at  the  two  discontinuities  are  nearly  equal  in  magnitude.  However, 
they  are  opposite  in  sign.  The  signal  waveform  seen  at  the  generator  end 
of  the  transducer  is  a superposition  of  these  two  reflections.  Let  us  express 
the  waveform  of  a step-function  voltage  which  ha?  travelled  a distance  2L 
down  a transmission  line  as  V(t)2L/  Then  the  signal  waveform  from  a 
small  disturbance  can  be  approximated  by 

Vg  > p[V(t)JL  - V(t-M)2L] 
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With  this  expression  we  can  readily  calculate  the  signal  waveform  if  we 
know  the  step-function  response  of  the  transmission  line. 

Figure  6 represents  the  signal  which  would  appear  at  the  input  of  a trans- 
ducer consisting  of  a 2.03  x 10  meter  radius  (No.  26)  copper  inner  con- 
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ductor  in  a 2.  41  x 10  meter  radius  copper  tube  (5  inch  refrigeration 
tubing).  There  is  a uniform  0.  5-millimeter  gap  between  the  inner  conduc- 
tor and  outer  conductor,  except  for  a one-meter  long  disturbed  section 
where  the  gap  is  0.49  millimeters,  located  1000  meters  from  the  driven 
end.  The  excitation  signal  is  a two-volt  step  function,  and  the  generator 
impedance  is  equal  to  the  nominal  cable  impedance.  Under  these  conditions, 
the  "TDR  picture"  would  be  a constant  1-volt  level  from  t = 0 until  t , the 
instant  where  the  first  reflection  from  the  disturbance  arrives  back  at  the 
generator  end.  For  an  air-insulated  line,  t would  be  6.  67  microseconds. 
Starting  at  t^,  the  waveform  of  Figure  6 would  be  added  to  the  constant 
1-volt  level.  (Note  that  the  polarity  of  this  superimposed  wave  is  negative 
because  the  disturbed  section  of  transducer  has  a lower  impedance  than  the 
rest  of  the  transmission  line.  ) This  waveform  corresponds  to  the  difference 
signal  which  appears  when  a disturbance  produces  a one-meter  long,  0.01- 
millimeter  displacement  in  the  position  of  the  inner  conductor. 

The  maximum  amplitude  of  the  signal  and  its  time  duration  are  shown  in 
Figure  6 as  \ anc*  T>  Actually,  the  total  time  duration  is  extremelv 
long,  but  most  of  the  energy  is  contained  between  the  half-voltage  points. 
Therefore,  we  have  used  this  definition  (full  width  at  half  maximum)  for  the 
time  duration  of  the  signal.  represents  the  effective  reflection 

coefficient;  t determines  the  time  resolution  and  thus  the  spatial  resolution 
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of  the  transducer.  For  this  example,  V,  is  about  15.6  microvolts 

MAX  _ 

(effective  reflection  coefficient  = 1.  56  x 10"^)  and  r is  6.2  x 10  seconds. 
The  corresponding  range  resolution  is  about  93  meters. 

To  calculate  the  waveform  of  a step-function  signal  which  has  propagated 
down  a long  transmission  line,  we  used  equations  developed  by  Wigington 
and  Nahman  (Reference  2).  These  equations  and  the  computation  procedures 
used  are  described  in  Appendix  C.  We  also  computed  signal  waveshapes  for 
a case  in  which  the  line  impedance  variation  is  expressed  as  a continuous 
analytical  function.  An  inverse  Fourier  transform  was  used  to  derive  the 
time  response  from  the  analytical  expressions  for  line  impedance  and  the 
transfer  function  of  a lossy  transmission  line.  Results  of  these  independent 
computations  are  included  in  Appendix  C.  They  show  that  the  step  approxi- 
mation for  the  conductor  displacement  is  indeed  valid  for  long  transducers, 
where  the  signal  waveform  is  determined  almost  entirely  by  the  transfer 
function  of  the  transmission  line.  On  a long  transducer,  the  amplitude  of 
the  reflected  signal  is  proportional  to  the  integral  of  the  displacement 
function  along  the  length  of  the  disturbance.  A 0.01  millimeter  displace- 
ment over  a 1 -meter  length  will  therefore  produce  almost  the  same  signal 
amplitude  and  duration  as  a 0.  001 -millimeter  displacement  over  a 10-meter 
length. 

Using  a double-step  approximation  for  the  inner  conductor  displacement,  we 
have  calculated  the  amplitude  and  time  widths  of  reflected  signals  in  uniform 
transducers  with  outer  conductor  ID's  of  0.  19,  0.  5 and  1 inch. 
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LENGTH  = 1000 


Figure  6.  Signal  Waveform;  0.19-Inch  Diameter,  1 km  Transducer 


To  show  the  effects  of  conductor  geometry  on  sensitivity  and  resolution,  we 

computed  and  plotted  V-T/.„  and  t versus  the  radius  ratio  (inner  conductor/ 

IVLnX 

outer  conductor)  for  various  conductor  gaps  and  for  outer  conductor  radii 
o 3 

of  2.4  x 10  meters  inch  refrigeration  tubing),  6.  35  x 10  meters  (£ 

_2 

inch  ID)  and  1.27  x 10  meters  (1  inch  ID).  In  all  cases,  the  disturbance 
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was  assumed  to  produce  a 0.  01 -millimeter  (3.94  x 10  inch)  deflection  of 
the  inner  conductor  over  a one-meter  section  of  the  transducer.  Lengths  of 
1000  meters  and  2000  meters  were  used  in  the  computations. 

.An  entire  family  of  curves,  similar  to  those  in  Figure  7,  has  been  gener- 
ated. In  Figure  7,  the  outer  conductor  diameter  is  \ inch,  the  gap  between 
inner  and  outer  conductor  is  0.05  millimeters,  and  the  disturbance  is  2000 
meters  from  the  driven  end. 

Optimum  Transducer  Design 

Both  the  effective  reflection  coefficient  and  the  time  duration  are  strongly 
influenced  by  the  line  attenuation.  In  all  of  the  cases  considered,  the  best 
sensitivity  and  resolution  occur  at  radius  ratios  of  approximately  0.3.  For 
larger  gaps,  where  the  inner  conductor  moves  closer  to  the  center,  the 
optimum  radius  ratio  tends  toward  0.  28,  which  is  the  minimum  attenuation 
condition  for  a concentric  transmission  line. 

Although  there  is  a clearly  defined  optimum  condition  for  radius  ratio,  the 
optimum  gap  between  inner  and  outer  conductor  is  not  so  easily  specified. 
There  is  a tradeoff  between  signal  amplitude  and  time  resolution  as  the 
gap  is  varied. 
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Figure  8 shows  the  effective  reflection  coefficient  and  range  resolution  of 
an  "optimum"  ^-inch  diameter  transducer  as  the  gap  is  varied  from  0.05 
to  5 millimeters.  For  a range  of  two  km.  the  gap  must  be  greater  than  one 
millimeter  to  achieve  a 30-meter  range  resolution.  Beyond  this  point,  the 
resolution  does  not  improve  very  rapidly  with  increasing  gap,  and  the  signal 
amplitude  continues  to  decrease  steadily. 

Predicted  Performance  versus  Length 

From  the  equations  describing  the  signal  waveforms,  we  can  readily  esti- 
mate the  effect  of  distance  on  the  signal  amplitude  and  time  duration.  At 
large  distances,  the  signal  amplitude  varies  inversely  as  the  square  of  the 
distance;  the  time  duration  increases  as  distance  squared.  This  relation- 
ship has  been  used  to  extrapolate  the  computed  points  and  generate  the 
curves  shown  in  Figures  9,  10  and  11.  These  curves  give  the  effective 
reflection  coefficient  and  range  resolution  versus  distance  for  0.  19,  0.  5 and 
1-inch  diameter  transducers  with  1-millimeter  nominal  gaps.  The  disturb- 
ance is  a 0.01  millimeter,  1 -meter  long  displacement  of  the  inner  conductor. 

If  we  choose  a range  resolution  of  30  meters  or  less,  and  an  effective  reflec- 
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tion  coefficient  of  3 x 10  or  more  as  the  acceptance  criteria  for  trans- 
ducers, we  can  use  Figures  9 through  11  to  establish  the  maximum  usable 
transducer  lengths.  A transducer  of  % -inch  refrigeration  tubing  with  a 
Number  16  inner  conductor  can  operate  satisfactorily  up  to  about  850  meters. 
With  an  outer  conductor  diameter  of  |inch,  lengths  up  to  two  km  are  possi- 
ble, and  a 1-inch  transducer  can  theoretically  be  deployed  in  lengths  up  to 
three  km. 


OUTER  DIAMETER  = 1/2  INCH 
OPTIMUM  INNER  CONDUCTOR  DIAMETER 
aG  = O.Ol  mm 
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Figure  8.  Reflection  Coefficient  and  Range  Resolution  versus  Gap 


Performance  Estimates  for  1-Inch  Transducer 


EFFECT  OF  NON-UNIFORM  TRANSDUCER  DIMENSIONS 


In  all  of  the  above  calculations,  a uniform  conductor  spacing  was  assumed. 

If  the  transducer  consists  of  a loosely  laid  wire  in  a tube,  the  spacing  will 
be  far  from  uniform.  This  will  increase  the  average  attenuation  and  will 
result  in  regions  of  low  sensitivity.  Furthermore,  multiple  reflections  of 
the  incident  step  waveform  and  the  return  signal  from  a disturbance  will 
result  from  the  impedance  variations  along  the  line  length.  This  causes 
additional  time  smearing  and  reduction  in  signal  amplitude.  It  also  results 
in  a form  of  crosstalk,  because  echoes  from  a disturbance  signal  will 
appear  in  all  subsequent  range  channels. 

Because  of  the  complexity  of  the  computations  involved,  it  was  not  practical 
to  consider  the  effects  of  multiple  reflections  and  frequency-dependent  losses 
at  the  same  time.  However,  we  did  perform  separate  computations  of  multi- 
ple reflections  on  a "lossless"  transmission  line  to  show  the  effect  of  dimen- 
sional variations  in  the  transducer. 

Electrical  Model  of  the  Loosely-Laid 
Wire-in-Tube  Transducer 


To  simulate  the  electrical  characteristics  of  the  loosely-laid  wire-in-tube 
sensor,  a computer-generated  curve  was  fitted  to  the  time-domain  reflec- 
tometer  measurements  on  the  10-meter  test  section  of  line.  We  considered 
deriving  a polynomial  expression  which  would  pass  through  selected  experi- 
mental data  points.  However,  this  degree  of  sophistication  was  not  neces- 
sary for  the  way  in  which  the  model  was  actually  used.  An  adequate  model 
was  produced  by  a linear  interpolation  method  in  which  the  data  points 
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entered  were  selected  to  provide  a reasonable  straight-line  fit  to  the  experi- 
mental curve.  Figure  12  is  a computer-generated  plot  of  the  reflection- 
coefficient  data.  The  inset  is  an  oscilloscope  camera  photograph  of  the  TDR 
display  (with  the  time  axis  reversed).  Figure  13  is  the  computed  impedance 
of  the  line,  based  upon  the  reflection-coefficient  data.  In  these  figures,  the 
horizontal  axis  is  shown  as  a time  scale.  A subsequent  measurement  was 
made  with  the  end  of  the  line  shorted  to  provide  a definite  marker  (multiple 
reflections  in  the  line  make  it  difficult  to  determine  the  exact  point  in 
Figures  12  and  13  where  the  10-meter  transducer  was  terminated  in  a 75- 
ohm  resistor).  It  was  found  that  the  round-trip  propagation  time  within  the 
10-meter  test  section  was  80  nanoseconds,  corresponding  to  a velocity 
factor  of  0.  833  with  respect  to  free-space  propagation.  To  relate  the  time 
scales  of  Figures  12  and  13  to  a physical  distance  along  the  transducer, 
note  that  zero  distance  on  the  transducer  corresponds  to  t = 5 nanoseconds, 
and  the  end  of  the  line  corresponds  to  t = 85  nanoseconds.  The  initial  5- 
nanosecond  delay  represents  a short  section  of  flexible  50-ohm  line  between 
the  TDR  and  transducer. 

Computation  Technique  for  Multiple  Reflections 

The  major  problem  of  generating  a model  for  the  lossless  non-uniform  wire 
transducer  is  to  develop  a mathematical  technique  for  incorporating  multi- 
ple impedance  variations  while  reducing  the  associated  bookkeeping  tasks. 
From  the  experimental  data  on  the  test  section  of  transducer  it  is  seen  that 
impedance  variations  occur  on  a spatial  scale  of  less  than  one  meter.  Anal- 
ysis of  a line  hundreds  of  meters  long  involves  hundreds  or  thousands  of 
discontinuities.  This  can  be  a staggering  bookkeeping  problem  if  the  mathe- 
matical method  is  poorly  chosen.  Memory  limitations  of  the  computer  further 
compound  the  problem. 
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TDR  DISPLAY 


By  treating  the  non-uniform  line  as  a series  of  short  line  segments  with 
equal  lengths,  the  time  history  of  the  propagated  and  reflected  wave  can  be 
kept  track  of  by  using  a method  proposed  by  Bewley  (Reference  3).  His 
method  is  graphically  represented  in  Figure  14.  Each  transmission  line 
junction  or  impedance  mismatch  can  be  expressed  as  a voltage  and  current 
level  at  that  junction.  The  equations  for  the  voltage  and  current  levels  for 
a given  junction  K are 


LK 


lk-l  ~ Vk+1  ~ Zk-1  !k-l  + Zk  Tk+1 
Zk+Zk-1 


and 


V.  = V.  . + Z.  , I.  , - Z.  ,1, 
k k-1  k-1  k-1  k-1  k 


The  above  equations  hold  for  all  junctions  other  than  the  special  cases  at 
the  generator  and  load  ends. 


The  equations  for  the  junction  between  the  generator  and  first  line  element 
are 
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impedance  of  the  generator 
source  voltage 


CIRCUIT 


Figure  14.  Bewley  Diagram 


The  junction  between  the  last  line  element  and  the  load  is  characterized  as 


and 


VL 


R1(VN  + ZNIN) 
RL  +ZN 


where 

R^  = load  impedance 
N = the  last  line  element 


Examination  of  the  previously  mentioned  equations  shows  the  dependence  of 
a given  junction  K on  its  left  and  right  junctions,  K-l  and  K+l  respectively. 
With  this  model,  keeping  a time  history  of  the  voltages  at  any  given  line 
junction  is  possible,  although  we  are  only  concerned  with  the  history  at  the 
generator  end. 

Two  more  observations  can  be  made  from  the  Bewley  diagram.  We  assume 
that  the  line  segments  are  of  equal  length,  and  that  the  velocity  of  propaga- 
tion of  the  traveling  wave  in  each  line  section  is  constant.  All  reflections 
occur  at  multiples  of  the  time  it  takes  to  traverse  one  line  segment.  We 
therefore  need  to  compute  reflections  only  at  these  discrete  intervals.  The 
second  observation  is  that  a traveling  wave  does  not  reach  junction  K until 
time  k.  Hence,  no  information  from  that  junction  is  received  at  the  gener- 
ator until  time  2k. 

Due  to  our  computer  limitations,  incorporation  of  the  math  model  into  a 
computer  package  was  a four  step  process.  The  four  steps  were: 
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1)  Generate  the  line  impedances 

2)  Plot  the  line  impedances  versus  line  element 

3)  Calculate  the  history  at  the  generator  end 

4)  Plot  the  history  versus  time 

Originally,  the  line  impedances  were  taken  from  the  experimental  time- 
domain  reflectometer  data  on  a 10-meter  transducer  section.  The  10-meter 
impedance  data  was  repeated  as  often  as  necessary  to  represent  the  length 
of  line  that  was  being  simulated.  This  scheme  turned  out  to  give  anomalous 
reflections  due  to  the  false  periodicity  induced  by  using  the  10-meter  TDR 
data.  In  actuality,  the  center  conductor  of  the  sensor  would  be  randomly 
lying  within  the  outer  conductor.  Therefore,  the  line  impedances  should  be 
generated  in  a random  manner  to  avoid  introducing  any  false  periodicity. 

To  avoid  generating  a line  that  would  have  large  changes  in  impedances 
from  one  segment  to  the  next,  a filter  was  devised.  This  filter  consisted  of 
averaging  four  consecutive  random  numbers.  This  eliminates  large  fluctu- 
ations in  impedance  values  between  adjacent  line  elements. 

Generation  of  the  N random  line  impedances  was  accomplished  in  the  follow- 
ing manner.  The  parameters  of  t.ie  transducer,  outer  conductor  radius, 
inner  conductor  radius,  and  inner  conductor  placement  with  respect  to  outer 
conductor  were  entered.  Filtered  random  numbers  were  drawn  to  determine 
the  plus  and  minus  offset  of  the  inner  conductor  with  respect  to  its  nominal 
placement  within  the  outer  conductor.  The  impedance  of  each  line  segment 
was  then  calculated  from  the  equation  for  an  eccentric  coaxial  transmission 
line.  Once  the  N line  impedances  were  generated,  they  were  written  onto 
magnetic  tape  for  later  use. 
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The  previously  described  algorithm  for  finding  the  time  history  at  the  gen- 
erator end  was  implemented  in  the  following  manner.  The  source  voltage, 
generator  impedance,  load  impedance,  and  line  impedances  which  were 
previously  calculated  and  stored  on  magnetic  tape  were  entered  into  the 
simulation  program.  It  should  be  noted  from  the  Bewlev  diagram.  Figure 
14,  that  for  each  time  interval  it  is  necessary  to  calculate  only  voltages 
and  currents  at  every  other  line  junction.  Hence,  at  t the  initial  voltage 
wave  is  calculated.  At  odd  numbered  times,  the  junctions  between  Elements 
1 and  2,  3 and  4,  etc.  are  updated,  and  at  even  times  the  junctions  between 
Elements  2 and  3,  4 and  5,  etc.  , are  updated.  The  junctions  at  the  gener- 
ator and  load  ends  are  updated  during  even  time  intervals. 

For  every  time  interval,  the  voltage  level  at  the  generator  and  the  corre- 
sponding time  interval  are  stored  in  a buffer.  After  200  time  intervals  have 
been  calculated,  the  buffer  is  then  written  on  tape.  The  process  is  then 
repeated  for  the  next  200  time  intervals  until  the  maximum  time  limit  is 
reached.  The  present  time  stop  is  set  to  2X  + 400  time  intervals.  This 
allows  the  wave  sufficient  time  to  travel  the  length  of  line  and  back  to  the 
generator. 

Since  the  algorithm  structure  requires  N line  impedances,  voltages  and  cur- 
rents to  be  retained  in  memory,  there  must  be  some  upper  limit  for  the 
value  of  N.  This  limit  will  vary  from  machine  to  machine,  so  in  our  case, 
using  a SDS  9300  computer,  the  value  for  N is  3200.  Without  restructuring 
the  model  or  using  peripheral  devices,  such  as  tape  drives  which  slow  down 
the  execution  of  the  program,  the  maximum  number  of  line  elements  that 
can  be  handled  is  3200. 


43 


In  addition  to  the  two  programs  that  were  just  described,  the  line  generator 
and  the  line  simulator,  there  are  two  more  programs  which  generate  a plot 
of  the  previously  calculated  data. 

The  first  routine  displays  the  generated  line  impedances,  stored  on  tape, 
versus  the  corresponding  line  element.  The  second  program  is  slightly 
more  complicated.  This  routine  has  two  options  available  to  the  user. 

The  first  option  is  to  read  a tape  containing  the  generator  history  informa- 
tion and  the  corresponding  time  interval,  then  graphically  display  that  infor- 
mation in  the  form  of  generator  history  versus  time.  The  second  option 
involves  reading  two  tapes  containing  two  different  sets  of  data,  comparing 
the  sets,  and  then  displaying  them  as  generator  history  versus  time.  The 
two  sets  of  data  are  the  results  of  two  line  simulation  runs  for  the  same 
length  transducer  where  one  run  has  one  or  more  perturbations  introduced 
to  the  line  for  simulation  of  intruders.  The  results  of  one  run  are  sub- 
tracted from  the  other,  thus  allowing  us  to  observe  the  time-varying  signals 
from  simulated  intrusions. 

Wire-in-Tube  Transducer  Simulation  Results 


Initial  computations  were  made  on  a simplified  model  of  the  transmission 
line  in  which  the  10-meter  section  was  broken  into  20  segments.  The  im- 
pedance of  each  segment  was  taken  as  the  average  value  from  the  experi- 
mental data  over  the  corresponding  interval.  A short-circuit  termination 
was  simulated  to  provide  information  on  the  fall  time.  It  was  hoped  that 
this  model  would  be  adequate,  because  the  solution  of  a 1-kilometer  line 
could  then  be  accomplished  with  a 2000  element  calculation  which  is  easily 


within  the  capabilities  of  our  computing  facility.  However,  the  numerical 
computation  of  the  reflected  voltage  for  this  model  showed  a fall  time  of 
zero  for  the  reflection  from  the  short-circuit  termination. 

Subsequent  computations  were  based  upon  a 200-element  approximation  to  the 
experimentally-determined  transducer  impedance.  The  computed  response 
showm  in  Figure  15  is  plotted  on  a normalized  time  scale,  where  each  time 
unit  represents  the  time  required  for  the  incident  wave  to  traverse  one  line 
segment.  Thus  it  takes  400  time  units  for  the  step-function  voltage  to 
travel  down  the  line  to  the  shorted  end  and  for  the  reflections  to  return  to  the 
input  end  of  the  200-segment  line. 

To  relate  this  to  the  TDR  data  shown  in  Figure  12,  each  time  unit  here 
corresponds  to  0.25  nanoseconds. 

The  computed  response  for  this  model  shows  that  the  fall  time  at  t = 400  is 
still  virtually  instantaneous.  In  other  respects,  the  response  appears  to  be 
correct  for  the  impedance  model  used. 

Neither  the  resistance  losses  nor  multiple  reflections,  as  computed  by  our 
numerical  technique,  could  explain  the  observed  fall  time.  Since  the  compu- 
tation technique  appeared  to  be  valid  in  other  respects,  and  the  200  element 
approximation  was  very  close  to  the  measured  line- impedance  variations  in 
the  test  transducer,  we  suspected  that  the  actual  measured  impedances  were 
incorrect.  We  know  that  the  inner  conductor  insulation  must  at  some  points 
be  touching  the  outer  tube,  and  at  other  points  the  inner  conductor  may  lie 
nearly  along  the  axis  of  the  outer  cylinder.  If  the  TDR  had  infinite  resolu- 
tion we  would  expect  the  impedance  to  go  from  a maximum  of  about  140  ohms 
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to  a minimum  of  50  or  60  ohms.  Since  the  impedances  derived  from  the 
TDR  measurements  only  ranged  from  about  70  to  105  ohms,  there  is  indeed 
a strong  reason  to  suspect  that  TDR  resolution  caused  some  "smoothing"  of 
the  measured  line  characteristics.  We  therefore  changed  the  model  by 
doubling  the  peak-to-peak  variations  about  the  85-ohm  nominal  impedance 
value.  The  computed  response  of  this  line  exhibited  a finite  fall  time,  but 
it  still  took  much  less  than  three  nanoseconds  to  reach  50  percent  of  the 
final  value. 

Taking  a further  critical  look  at  the  measured  impedance  data,  we  noted  that 
the  far  end  of  the  transducer  seemed  to  have  fewer  and  less  extreme  imped- 
ance variations  than  the  first  few  meters  of  line.  There  is  no  reason  why 
this  should  be  true;  in  fact  the  inner  conductor  must  touch  the  wall  at  least 
every  meter  or  so  along  the  entire  line  because  it  was  not  under  tension. 
Evidently  the  waveform  degradation  of  the  line  itself  was  introducing  addi- 
tional smoothing  after  the  first  few  meters  of  line. 

To  correct  for  this  effect,  we  simply  took  the  impedance  data  from  the  first 
few  meters  of  transducer  and  repeated  it  periodically  until  we  had  "extended" 
this  new  model  to  200  segments.  Finally,  we  doubled  the  impedance  vari- 
ations to  produce  the  model  shown  in  Figure  16;  the  computed  response  of 
this  line  is  given  in  Figure  17.  Nowr  the  50  percent  fall  time  lies  between 
10  and  20  units  (2.  5 to  5 nanoseconds)  and  the  reflected  waveform  shows  the 
same  smoothing  effect  on  the  last  half  of  the  line  as  we  see  in  the  TDR  data. 
This  final  model  is  probably  a fairly  accurate  simulation  of  the  real  imped- 
ance variations  in  a loosely-laid  wire  transducer,  and  the  w:aveform  degrad- 
ation in  such  a transducer  is  predicted  by  the  numerical  computation. 
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We  extended  the  above  model  to  1600  elements  and  calculated  its  step 
response.  When  the  generator  impedance  is  matched  to  the  average  trans- 
ducer impedance  of  about  85  ohms,  and  the  transducer  is  terminated  with  an 
85  ohm  resistor,  the  simulated  "TDR  picture"  of  a 1600-element  line  is 
shown  in  the  six  successive  simulation  plots  of  Figure  18.  A substantial 
loss  in  amplitude  of  reflections  is  evident  for  the  far  end  of  the  line.  (Note 
the  vertical  scale  changes  on  the  various  portions  of  the  line.  ) The  trans- 
ducer model  consists  of  a basic  50  segment  line,  repeated  32  times.  There 
appear  to  be  anomalous  responses  at  about  t = 1400  and  t = 2800.  These 
are  probably  artifacts  of  the  periodic  line  formed  by  repeating  the  imped- 
ance variations  every  50  segments.  Multiple  reflections  can  add  construc- 
tively or  destructively  to  produce  reinforcement  or  cancellation  at  certain 
points  along  the  line. 

To  simulate  the  effects  of  time-varying  disturbances,  we  computed  the 
response  of  an  undisturbed  line,  then  introduced  impedance  variations  in 
the  model  at  selected  points  and  recomputed  the  response,  and  finally  plot- 
ted the  difference  between  the  two  runs.  The  disturbances  were  simulated 
by  adding  1.  7 ohms  to  the  initial  values  of  line  impedance  over  selected 
20-segment  intervals.  This  corresponds  to  a disturbance  which  produces 
a nominal  reflection  coefficient  change  of  0.01  over  a 1 -meter  interval. 

The  six  successive  simulation  plots  of  Figure  19  show  the  results  of  a 
multiple-disturbance  run.  The  first  disturbance,  which  begins  at  segment 
number  10  (t  = 20),  shows  a distinct  response  with  an  average  value  of 
about  0.01.  Following  this  disturbance,  we  see  the  "echoes"  from  subse- 
quent impedance  variations  in  the  line.  The  amplitude  of  the  response  from 
other  simultaneous  disturbances  beginning  at  segments  number  110  and  310 
(t  = 220  and  620)  is  not  appreciably  diminished,  but  there  is  a noticeable 
broadening  of  the  time  response. 
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T igure  18.  Simulated  TDK  Picture  for  1600-Segment  Loosely-Paid  Wire  Transducer 
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Loosely-Laid  Wire  Transducer  (continued) 


85  OHMS  SOURCE  VOLTAGE  = 2.0  VOLTS 


Figure  18.  Simulated  TDR  Picture  for  1600-Segment 
Loosely-Laid  Wire  Transducer  (continued) 
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Figure  18.  Simulated  TOR  Picture  for  1600-Segment 

Loosely-Laid  Wire  Transducer  (concluded) 
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Transducer  with  Multiple  Reflections  (continued) 


Varying  Response  of  Loosely-Laid  Wire 
Multiple  Reflections  (continued) 
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Transducer  with  Multiple  Reflections  (concluded) 


Again,  the  line  shows  anomalous  behavior  at  about  t = 1400.  The  calculated 
response  shows  time- varying  reflections  from  undisturbed  portions  of  the 
line  which  are  comparable  in  magnitude  to  the  signals  from  actual  distur- 
bances. The  last  disturbance  at  segment  number  1510  (t  = 3020)  is  difficult 
to  discern  among  the  echoes  from  the  other  disturbances. 

In  an  attempt  to  remove  the  anomalous  behavior  at  t = 1400  and  t = 2800,  we 
tried  changing  the  model  to  increase  the  fundamental  length  from  50  to  200 
segments.  This  was  done  by  using  a time  inversion  of  the  first  50  segments 
between  50  and  100,  an  inversion  about  the  average  value  of  impedance  for 
segments  100  through  150,  and  another  time  inversion  to  generate  segments 
150-200.  Although  the  resulting  response  was  altered,  there  were  still 
regions  of  anomalous  reflections. 

We  did  not  make  further  refinements  of  the  loosely-laid,  enclosed-wire 
transducer  model  because  this  transducer  appears  far  from  optimum  for 
pulsed-mode  operation.  The  simulated  1600-element  transducer  repre- 
sents a physical  length  of  80  meters.  Even  at  this  length,  attenuation  and 
crosstalk  problems  are  severe.  Without  a highly  sophisticated  processor, 
the  maximum  usable  length  of  a loosely-laid  wire-in-tube  transducer  will 
probably  be  less  than  a few'  hundred  meters. 

Effect  of  Random  Gap  Variations  in 
a 0.  5-Inch  Optimized  Transducer 

Further  computations  of  multiple- reflection  effects  were  carried  out  for  a 
0.  5-inch  diameter  transducer.  The  inner  conductor  diameter  was  chosen 
as  0.  14  inch  and  the  nominal  gap  was  set  at  one  millimeter.  These 
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dimensions  were  shown  earlier  to  provide  nearly  optimum  electrical  charac- 
teristics for  a uniform  transducer.  To  avoid  the  anomalous  results  obtained 
with  periodic  impedance  variations,  we  used  a transducer  model  in  which 
the  dimensional  variations  are  completely  random.  This  model  was  created 
by  using  a random-number  generator  routine  in  the  computer  to  represent 
the  inner  conductor  displacement.  The  inner  conductor  displacement  cannot 
vary  abruptly  between  its  extreme  positions  in  a real  transducer.  We  there- 
fore used  a "filtering"  process  consisting  simply  of  a four-number  running 
average  to  smooth  the  random-number  generator  output. 

For  the  first  simulation  we  used  a 1600-element  transducer  with  peak-to- 
peak  fluctuations  of  about  0.  7 millimeters  about  the  1 -millimeter  nominal 
gap.  Figure  20  is  a plot  of  the  first  300  elements  in  this  transducer  model 
showing  the  variation  in  gap. 

To  simulate  the  effect  of  a time-varying  disturbance,  we  first  calculated 
the  time  history  of  the  voltage  at  the  generator  end  of  the  transducer.  Then 
the  transducer  model  was  altered  by  adding  0.  01  millimeter  to  the  gap  for 
Elements  100  through  103,  and  Elements  1500  through  1503.  The  voltage 
response  was  computed  again,  and  the  "undisturbed"  response  then  sub- 
tracted. When  plotted,  the  result  is  a time-amplitude  history  showing,  in 
effect,  the  outputs  of  a very  large  number  of  processed  range  channels  when 
two  similar  time- varying  disturbances  are  present  simultaneously. 

Figure  21  (two  parts)  shows  portions  of  the  time-amplitude  plot  for  the 
1600-element  transducer.  The  round-trip  time  delay  for  each  transducer 
element  is  two  time  units.  We  can  think  of  each  time  unit  as  a separate 
time-gated  range  channel,  with  the  corresponding  range  increment  equal 
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Figure  21.  Partial  Time-Amplitude  History  for 
1600-Element  Transducer  (concluded) 
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to  one-half  the  length  of  a transducer  element.  We  will  discuss  the  physi- 
cal significance  of  these  arbitrary  time-unit  scales  later.  The  response 
from  the  disturbance  at  100  elements  appears  at  t = 200,  and  the  disturbance 
at  1500  elements  appears  at  t = 3000. 

This  transducer  does  not  exhibit  a serious  variation  in  sensitivity  or  degrad- 
ation in  time  response  along  its  length.  Echoes  of  the  disturbances  are 
seen  in  subsequent  range  channels,  but  their  levels  are  relatively  low.  In 
this  simulation,  the  transducer  was  terminated  in  a short  circuit.  The 
position  of  the  short  is  not  distinctly  evident  at  t = 3200,  but  there  is  a 
strong  echo  200  time  units  later,  at  t = 3400.  This  is  a reflection  of  the 
disturbance  signal  which  appeared  at  t = 200,  after  it  has  travelled  down 
the  line  and  reflected  from  the  short.  Its  amplitude  has  decreased  by  a 
factor  of  two  and  it  is  somewhat  broadened  in  time.  Time  units  beyond 
t = 3200  are  only  of  academic  interest  because  all  of  the  primary  signals 
have  already  been  observed.  It  is  obvious,  however,  that  the  transducer 
should  be  terminated  with  a matched  load  so  that  persisting  echoes  will  be 
as  small  as  possible  when  the  next  step  interrogates  the  transducer. 

A further  check  on  the  effects  of  dimensional  variations  was  made  by  com- 
puting the  output  of  a similar  ^-inch  diameter  transducer.  The  same 
"random"  function  was  used  to  describe  the  gap  variations  as  in  the  case 
just  described.  However,  the  magnitude  of  the  variations  was  doubled,  so 
that  the  peak-to-peak  deviations  in  the  inner-conductor  position  are  approxi- 
mately 1.  5 millimeters.  The  transducer  length  was  also  extended  to  3200 
elements.  Disturbances  were  introduced  at  100,  1500  and  3000  elements. 
Using  the  same  computations  as  before.  Figure  22  (three  parts)  shows 
three  time-amplitude  plots  in  the  region  of  interest. 
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T igure  22.  Partial  Time-Amplitude  History  for  3200-Element  Transducer 
with  1. 5 mm  Gap  Variations  (continued) 
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The  first  disturbance  signal,  at  t = 200,  is  very  similar  to  the  previous 
case.  However,  the  echoes  following  the  disturbance  are  larger  in  ampli- 
tude by  approximately  a factor  of  four.  A substantial  broadening  in  time 
and  reduction  in  amplitude  is  evident  for  the  second  disturbance  signal  at 
t = 3000.  There  is  also  an  apparent  delay  before  the  signal  reaches  its 
maximum  amplitude.  All  of  these  effects  are  even  more  pronounced  at 
t = 6000,  where  the  signal  no  longer  stands  out  clearly  from  the  "noise". 

The  noise  is  actually  caused  by  multiple  echoes  of  the  other  signals. 

When  presented  in  the  above  form,  the  results  are  overly  pessimistic. 

They  represent  an  extremely  high-resolution  picture  of  the  time  history, 
which  would  not  be  practical  or  even  useful.  A real-world  signal  processor 
will  have  only  a limited  number  of  range  channels,  each  of  which  forms  an 
average  over  many  time  units.  This  will  substantially  smooth  the  "noise,  " 
which  has  a long-term  average  value  of  zero.  The  effect  of  this  smoothing 
has  been  determined  by  computing  a coarse-resolution  version  of  the  time 
amplitude  history  that  was  partially  depicted  in  the  three  plots  of  Figure  22. 

- 

Signal  amplitudes  were  averaged  over  25  consecutive  intervals  of  250  time 
units.  The  result,  shown  in  Figure  23,  simulates  the  outputs  of  a 25-channel 
processor  with  range  gate  widths  of  250  time  units.  The  three  channels 
\ containing  actual  disturbance  signals  stand  out  clearly  above  the  worst-case 

l,  crosstalk  level  (note  that  the  bottom  level  of  the  histogram  is  actually  a 

negative  voltage).  However,  it  is  still  apparent  that  excessive  dimensional 
variations  in  the  transducer  will  limit  the  dynamic  range  of  the  sensor  be- 
cause  of  echoes  or  crosstalk.  The  crosstalk  can  be  suppressed  by  signal 
processing,  but  a relatively  sophisticated  processor  might  be  required. 
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Figure  23.  Simulated  Processor  Outputs  for  3200-Element  Transduc 


All  of  the  multiple- reflection  computations  for  the  0.  5-inch  transducer  have 
been  presented  on  an  arbitrary  scale  of  time  units.  This  is  because,  for  a 
lossless  transducer,  the  actual  physical  length  (and  therefore  the  time  de- 
lay) has  no  influence  on  the  waveform.  The  only  parameters  of  importance 
are  the  number  of  impedance  variations  encountered  and  their  magnitudes. 
To  relate  these  results  to  actual  transducer  lengths,  we  must  define  the 
distance  scale  over  which  the  dimension  variations  occur. 

In  the  case  of  the  loosely-laid  #26  wire  in  a 4 -inch  tube,  the  TDR  measure- 
ments showed  that  large  impedance  variations  took  place  over  short  dis- 
tances. On  the  average,  there  were  several  "cycles"  of  impedance  vari- 
ation per  meter.  With  larger  conductors  which  have  less  tendency  to 
deform  over  short  distances,  the  spatial  period  of  impedance  variations  can 
possibly  be  increased  to  an  average  of  one  or  two  meters.  Referring  to 
Figure  20,  there  are  about  65  "cycles"  of  impedance  variation  in  the  first 
400  elements  of  the  transducer  model.  If  an  average  cycle  corresponds  to 
1.5  meters,  the  equivalent  physical  length  represented  by  the  400-element 
picture  is  about  100  meters.  Each  meter  of  transducer  length  then  corre- 
sponds to  four  transducer  elements  or  eight  time  units.  With  this  assump- 
tion concerning  the  spatial  variations  in  transducer  impedance,  the  3200- 
element  case  would  represent  an  800-meter  transducer.  Disturbances  are 
located  at  25,  375  and  775  meters.  In  Figure  23,  each  range  channel  is  31 
meters  wide. 


SIGNAL-TO-NOISE  ESTIMATES 


We  have  made  preliminary  estimates  of  the  signal-to-noise  ratios  attainable 
with  the  pulsed-mode,  enclosed-wire  transducer.  These  estimates  are 
based  upon  the  simulation  results  for  a uniform  transducer,  and  show  the 
factors  limiting  the  signal-to-noise  ratio.  The  calculations  show  that  even 
a simple  processor  can  give  an  acceptable  ratio  of  signal  to  the  theoretical 
noise  floor  of  the  electronics. 

For  the  signal-to-noise  calculations,  we  assumed  a single-channel  processor 

of  the  form  showm  in  Figure  24.  It  consists  simply  of  a pulse  generator 

which  periodically  applies  a step-function  drive  signal  to  the  transducer. 

The  generator  impedance  is  matched  to  the  transducer  impedance  Zq 

to  minimize  reflections.  A video  preamplifier  is  incorporated  to  overcome 

the  noise  contributions  of  subsequent  switching  and  amplifier  circuits.  To 

simplify  the  noise  calculations,  we  will  assume  that  the  preamplifier  has 

unity  gain,  and  that  the  dominant  noise  sources  are  the  thermal  noise  in  the 

transducer  and  the  equivalent  input  noise  of  the  preamplifier.  Following 

the  video  preamplifier  is  a range-gate  switch  and  an  integrating  capacitor. 

The  range-gate  switch  closes  for  a time  interval  At  after  a delay  time  t 

following  the  drive  pulse.  If  the  velocity  of  propagation  in  the  transducer 

is  V,  the  range  gate  selects  the  reflected  signals  from  the  range  interval 

AR  = VAt/2  located  at  a distance  X = Vt  / 2. 

o 

The  integrating  capacitor  takes  the  average  value  of  signal  plus  noise 
voltage  during  the  time  the  range-gate  switch  is  closed,  and  holds  this 
value  while  the  switch  is  open.  A bandpass  amplifier  following  the  range- 
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gate  circuit  limits  the  processor  bandwidth  to  the  desired  signal  bandwidth 
B.  We  will  assume  that  the  range  gate  interval  is  centered  on  the  signal 
pulse. 

If  the  drive  voltage  is  E and  the  effective  reflection  coefficient  of  the  dis- 
turbed transducer  section  is  o^pp,  the  signal  voltage  has  a peak  value  equal 
to  Ep,,™.  The  exact  value  of  signal  voltages  can  be  obtained  by  integrating 
the  signal  pulse  over  the  time  interval  At.  When  the  range-gate  interval  is 
centered  on  the  signal  pulse,  we  can  use  simplified  expressions  which  give 
fairly  accurate  approximations  for  the  average  signal  voltage.  If  the  half- 
voltage width  t of  the  signal  pulse  is  greater  than  At,  the  average  value  of 
signal  voltage  is  close  to  the  peak  value 
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Figure  24.  Processor  Assumed  in  Signal-to-Noise  Calculations 
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When  the  range-gate  interval  is  longer  than  the  pulse  width,  the  average 
signal  voltage  is  reduced  by  the  factor  K^.  The  value  of  K is  approxi- 
mately unity,  and  e^  can  be  approximated  by 

6s  - E ^FF  At  (T<  At) 


The  rms  noise  voltage  e from  a source  at  temperature  T and 


resistance 


Zq  in  a bandwidth  Af  is 


6no  = k T0  ZQ  Af)?  volts  rms 


where  k TQ  = 4.  14  x 10-“  joules  at  Tq  = 300°K.  A low-noise  video  pre- 
amplifier will  add  less  than  three  dB  to  the  thermal  noise  limit.  The  effec- 
tive noise  bandwidth  is  approximately  1/At.  Assuming  a 3-dB  noise  figure 
for  the  video  preamplifier,  we  have 


/8  k T ZV 

6n  =v  ir-^)  volts  rms 


Using  radar  terminology,  the  "single-pulse"  signal-to-noise  voltage  ratio 
is  approximately 
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These  expressions  give  the  ratio  of  peak  signal  to  rms  noise  voltage  at  the 
low-frequency  amplifier  input.  The  best  signal-io-noise  ratio  is  obtained 
when  At  is  "matched"  to  the  pulse  width  t. 


An  improvement  in  output  signal-to-noise  ratio  is  gained  if  the  bandpass 
amplifier  has  a bandwidth  less  than  the  maximum  bandwidth  allowed  by  the 
Nyquist  sampling  criterion 

B -2SH 
max  2 


where  PRF  is  the  pulse  repetition  frequency.  The  signal  voltage  is  propor- 
tional to  B /B,  and  the  noise  voltage  is  proportional  to  (B  )^/B.  Thus 
max  max  « 

there  is  an  "integration  gain"  in  the  signal-to-noise  ratio  equal  to  (PRF/2B)  . 


The  maximum  pulse  repetition  frequency  is  determined  by  the  total  trans- 
ducer length  L and  the  velocity  of  propagation  V. 


PRF 


max 


V 

2L 


We  can  also  express  the  optimum  gate  interval  At  = t in  terms  of  the  trans- 
ducer range  resolution  AR 

T = 2AR 


Combining  these  factors,  the  output  signal-to-noise  voltage  ratio  when  the 
processor  gate  time  is  optimized  is  approximately 
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Expressed  in  dB,  the  optimum  signal-to-noise  ratio  is 

<5  (E  AR 

o - i n 1 EFF 

N " Og10  1 6kT  Z BL 
opt  o o 

A numerical  example  of  the  realizable  signal-to-noise  ratio  with  an  opti- 
mized transducer  and  processor  is  given  here  for  the  case  of  a 2 km,  \- 
inch  transducer  with  a uniform  1 -millimeter  gap.  Assume  a drive  signal 
E = 1 volt  and  the  signal  bandwddth  B = 1 kHz.  The  transducer  impedance 
Zq  is  approximately  50  ohms.  From  the  transducer  performance  curves, 
we  find  that  at  2 km,  AR  = 30  meters  and  = 3 x 10"°.  Under  these 

i^r  r 

conditions,  the  output  signal-to-noise  ratio  is  approximately  36  dB. 

The  above  signal-to-noise  ratio  expressions  relate  the  peak  signal  level  to 
the  rms  noise  voltage.  In  an  intrusion  sensor  with  a threshold-detector 
alarm  circuit,  we  are  concerned  with  the  peak  values  of  both  the  noise  and 
signal  voltage.  For  Gaussian  noise,  the  peak  value  will  very  rarely  exceed 
three  times  the  rms  value.  The  signal-to-noise  voltage  ratios  given  above 
can  be  divided  by  three  to  give  a more  conservative  estimate  for  false- 
alarm  rate  predictions.  This  reduces  the  signal-to-noise  ratio  in  the  above 
example  to  26  dB,  which  is  still  adequate  to  provide  a very  low  false-alarm 
rate. 
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In  the  calculations  presented  above,  we  have  not  taken  into  account  the 
cumulative  effects  of  background  noise,  fluctuations  in  the  drive  signal 
level,  or  the  "clutter"  effects  caused  by  multiple  reflections.  The  equations 
indicate  that  the  signal-to-noise  ratio  can  be  improved  indefinitely,  simply 
by  increasing  the  drive  signal  level.  However,  the  background  noise  due  to 

I actual  seismic/acoustic  disturbances  increases  in  direct  proportion  to  the 

drive  level.  Fluctuations  in  the  drive  signal  level  can  easily  exceed  a few 
parts  in  10®,  which  sets  a lower  limit  of  about  10  ® on  the  minimum  effec- 
tive reflection  coefficient  for  reliable  detection. 

In  the  simple  processor  circuit  assumed  for  the  signal-to-noise  calculations, 
the  dynamic  range  of  the  video  preamplifier  and  range-gating  circuits  sets 
a limit  on  the  allowable  drive  level.  This  difficulty  can  be  avoided  with  a 
slight  increase  in  circuit  complexity,  as  discussed  later. 

Clutter  due  to  multiple  reflections  in  a non-uniform  transducer  is  a station- 
ary effect,  and  is  not  inherently  a source  of  noise.  However,  crosstalk 
from  multiple  reflections  can  cause  background  noise  along  the  entire 
transducer  length  to  contribute  to  the  noise  level  in  each  range  channel. 

We  expect  the  cumulative  background  noise  voltage  from  incoherent  sources 
to  be  proportional  to  the  square  root  of  the  number  of  range  channels,  multi- 
plied by  the  average  interchannel  crosstalk  level.  As  long  as  the  crosstalk 
level  is  low,  the  cumulative  noise  contribution  will  not  be  severe.  With  a 
3-km  transducer  and  100  30-meter  range  channels,  the  cumulative  back- 
ground noise  due  to  crosstalk  from  all  other  channels  will  not  exceed  the 
background  noise  in  a given  channel  if  the  crosstalk  level  is  below  20  dB. 
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Unless  the  transducer  is  extremely  uniform,  the  signal-to-clutter  ratios 
will  be  low.  For  example,  if  the  displacement  caused  by  an  intruder  is 
0.  01  millimeter,  and  there  is  a 0.  3 millimeter  variation  in  the  nominal 
gap,  the  signal-to-clutter  ratio  will  be  on  the  order  of  -30  dB. 

Excitation  Waveforms 

Our  sensor  analysis  was  based  upon  the  assumption  of  a fast-rise  step- 
function  drive  signal.  This  is  the  basic  excitation  signal  used  in  time- 
domain  reflectometry,  and  it  is  convenient  for  time-domain  as  well  as 
frequency-domain  analysis.  However,  it  is  not  necessarily  the  optimum 
excitation  for  the  pulsed-mode  transducer.  We  did  not  perform  an  analysis 
of  more  complex  excitation  signals,  but  we  identified  variations  of  the  step 
function  which  offer  performance  advantages. 

In  the  simple  case  of  a step-function  drive  signal,  we  must  bring  the  drive 
voltage  back  to  zero  before  another  step  can  be  applied.  Reflections  from 
the  transducer  will  also  occur  from  this  "retrace"  signal.  Rather  than 
accepting  the  additional  50  percent  dead  time  required  for  these  reflections 
to  damp  out,  we  can  just  as  well  use  a dual-slope,  step-function  excitation- - 
in  other  words,  a square  wave.  All  that  is  necessary  is  to  reverse  the 
polarity  of  the  return  signals  during  alternate  half-cycles  of  the  excitation 
square  wave  so  that  the  return  signals  will  be  combined  constructively. 

When  square-wave  excitation  is  used,  the  actual  frequency  of  the  drive 
signal  is  one-half  the  effective  PRF. 
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This  square-wave  form  of  excitation  has  another  benefit  in  addition  to 
reduction  of  dead  time.  The  drive  signal  can  be  made  bipolar,  with  an 
average  value  of  zero.  This  removes  any  unintentional  d-c  bias  from  the 
transducer.  Since  the  polarity  of  the  detection  electronics  is  alternated  at 
a high  rate,  extraneous  d-c  or  low-frequency  signals  are  cancelled.  These 
features  make  the  pulsed-mode  transducer  insensitive  to  the  basic  low- 
frequency  signals  generated  by  the  wire-in- tube  transducer  when  it  is 
biased,  either  externally  or  by  electret  effects  in  the  insulation. 

A simple  implementation  of  the  square-wave  drive  technique  is  shown 
functionally  in  Figure  25.  The  only  additions  to  the  basic  circuit  are  a 
differential  video  preamplifier,  a reference  load  resistance,  and  four 
switching  diodes. 

Synchronous  polarity  switching  in  the  detection  electronics  is  automatically 
provided  by  a diode  ring  in  combination  with  the  bipolar  excitation  signal. 
When  the  drive  signal  polarity  is  positive,  the  transducer  is  connected  to 
the  non-inverting  input  of  the  differential  video  amplifier  and  the  reference 
load  is  connected  to  the  inverting  input.  A negative  drive  signal  reverses 
these  connections. 

As  shown  in  the  diagram,  the  generator  and  reference  load  impedances  are 
equal  to  the  nominal  transducer  impedance.  It  is  necessary  to  match  the 
generator  impedance  to  the  transducer  to  suppress  multiple  reflections. 

A reference  load  with  the  same  impedance  as  the  transducer  provides 
suppression  of  noise  which  may  be  present  in  the  drive  signal  itself.  It  also 
reduces  the  dynamic-range  requirement  of  the  video  preamplifier  and  sub- 
sequent range-gating  electronics,  which  would  otherwise  have  to  accommo- 
date the  full  drive-signal  amplitude  multiplied  by  the  preamplifier  gain. 
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Figure  25.  Square-Wave  Drive  Technique 


There  is  a slight  penalty  in  signal-to-noise  ratio  with  this  circuit  because 
of  the  additional  amplifier  input  and  the  shot-noise  contribution  of  the 
switching  diodes.  Calculations  show  that  the  combined  effect  can  be  less 
than  6 dB,  which  can  be  made  up  by  doubling  the  drive  voltage.  Substan- 
tially higher  drive  voltages  can  be  used  with  this  differential  system  than 
with  a single-ended  system,  as  long  as  the  common-mode  input  limits  of 
the  differential  preamplifier  are  not  exceeded. 

The  risetime  of  the  excitation  signal  was  assumed  to  be  zero  in  our  trans- 
ducer analyses.  A fast-risetime  signal  is  essential  for  high-resolution 
time-domain  reflectometry.  However,  the  resolution  of  the  pulsed-mode 
enclosed-wire  sensor  is  limited  by  waveform  degradation  in  the  transducer, 
as  well  as  by  practical  considerations  in  the  processor.  Using  an  excita- 
tion signal  with  a risetime  much  shorter  than  the  resolution  limits  of  the 
transducer  and  processor  does  not  offer  any  real  advantage.  If  we  increase 
the  drive-signal  risetime,  the  peak  amplitudes  of  short-range  signal  and 
clutter  reflections  at  the  preamplifier  input  will  be  reduced.  However,  the 
signal  amplitude  at  the  processor  output  will  not  be  affected  until  the  rise- 
time of  the  drive  signal  approaches  the  range  gate  time  interval.  An  exci- 
tation signal  with  a risetime  on  the  order  of  one-half  the  gate-time  interval 
will  reduce  the  dynamic-range  requirements  of  the  video  preamplifier  with- 
out appreciable  degradation  in  sensor  resolution. 

Alternate  Signal- Process  or  Designs 

Many  of  the  waveforms  and  processing  techniques  used  in  modern  radar  and 
sonar  systems  may  be  applicable  to  the  pulsed-mode  enclosed-wire  sensor. 
The  signal-to-noise  estimates  show  that  a simple  processor  can  be  used 


with  transducer  lengths  of  two  km  or  more.  Because  of  the  presence  of 
clutter,  crosstalk  and  background  noise,  it  is  difficult  to  say  how  much 
improvement  could  be  gained  with  a more  complex  driving  waveform  and  a 
sophisticated  signal  processor.  Existing  analyses  of  the  pulsed-mode  trans- 
ducer need  to  be  supplemented  with  experimental  data  and  proof  of  feasibility 
before  the  detailed  analysis  and  design  of  a complex  processor  is  justified. 

One  type  of  processor  which  could  offer  a significant  cost  advantage  for  a 
multichannel  sensor  is  based  upon  charge-coupled  devices  (CCDs).  The 
functions  of  a k-channel  processor  can  be  implemented  with  time-shared 
recursive  filters  using  k-bit  CCD  analog  delay  lines.  Outputs  from  all 
channels  can  be  time-multiplexed  to  appear  on  a single  line.  This  would 
eliminate  the  need  for  separate  range  switches  and  analog  filters  in  each 
channel. 

SUGGESTED  TRANSDUCER  CONFIGURATIONS 

Our  analysis  shows  that  the  loosely-laid  wire-in-tube  configuration  is  not 
well  suited  for  pulsed-mode  operation.  Since  the  inner  conductor  can 
assume  any  position  within  the  outer  tube,  the  transducer  impedance  vari- 
ations are  very  large.  Resulting  multiple  reflections  cause  severe  attenu- 
ation, clutter  and  crosstalk.  An  alternate  transducer  design  must  be  used 
unless  the  maximum  length  is  less  than  a few  hundred  meters.  It  is  diffi- 
cult to  satisfy  all  of  the  conflicting  design  requirements.  The  inner  conduc- 
tor must  be  spaced  fairly  uniformly  within  the  outer  tube.  Peak-to-peak 
variations  in  the  gap  should  not  exceed  approximately  1 millimeter.  But  in 
order  to  preserve  the  sensitivity  of  the  wire-in-tube  transducer,  the  inner 
conductor  must  be  able  to  move  freely.  It  is  desirable  to  preserve  the 
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mechanical  resonances  which  are  responsible  for  high  sensitivity  in  the 
"strumming"  mode.  The  pulsed-mode  transducer  should  also  be  cheap, 
rugged  and  easily  deployed.  Transducer  configurations  which  show  promise 
of  meeting  these  requirements  are  identified  in  this  section. 

A good  compromise  between  cost  and  performance  can  be  achieved  with  a 
£-inch  diameter  aluminum  outer  conductor.  The  analysis  shows  that  a 
signal-to-noise  ratio  of  25  dB  and  a range  resolution  of  30  meters  is  theo- 
retically possible  with  a 2-km  long,  |-inch  diameter  transducer.  One  inch 
tubing  provides  comparable  performance  at  50  percent  greater  distances, 
but  costs  and  weighs  about  three  times  as  much  as  ^-inch  tubing.  Alumi- 
num has  a slightly  higher  resistivity  than  copper,  but  for  high-conductivity 
alloys  the  difference  is  only  about  30  percent.  High-frequency  resistance 
is  proportional  to  the  square  root  of  resistivity  and  inversely  proportional 
to  conductor  diameter.  Resistive  losses  in  the  transducer  will  be  domin- 
ated by  the  inner  conductor,  and  the  choice  of  copper  or  aluminum  for  the 
outer  conductor  will  affect  the  total  resistance  by  only  about  4 percent. 

Aluminum  tubing  also  costs  much  less  than  copper.  Half-inch  copper 
refrigeration  tubing  costs  about  50  cents  per  foot;  aluminum  tubing  is  less 
than  15  cents  per  foot.  Another  advantage  of  aluminum  is  its  weight.  Five 
hundred  meters  of  ^-inch  aluminum  tubing  weighs  less  than  100  pounds; 
copper  tubing  would  weigh  nearly  300  pounds.  Fabrication  may  be  some- 
what more  difficult  with  aluminum  than  copper,  although  weather-proof 
joints  and  connectors  are  commonly  used  in  aluminum-jacketed  coaxial 
lines  for  cable  TV  systems. 


It  is  also  possible  to  use  a braided  copper  outer  conductor  with  a supporting 
dielectric.  This  may  offer  advantages  in  fabrication  and  handling. 

The  most  critical  consideration  in  the  transducer  design  is  the  inner  con- 
ductor and  the  way  in  which  it  is  supported.  For  a \ -inch  outer  conductor 
the  optimum  inner  conductor  diameter  is  0.  14  inches  (equivalent  to  #7  wire). 

The  inner  conductor  diameter  can  be  reduced  to  0.  075  inches  before  the 
attenuation  becomes  excessive.  This  is  approximately  the  diameter  of  a 
#12  wire.  We  list  below  a few  of  the  options  which  can  be  considered  for 
the  inner  conductor  and  its  support  structure. 

Stranded  Wires--Solid  #12  or  larger  wire  is  too  rigid  to  lie  in  the 
outer  tube  without  extreme  variations  in  gap.  However,  multiple- 
strand  wire  may  be  sufficiently  flexible  to  conform  to  the  contour 
of  a deployed  £ -inch  tube. 

Braided  Conductors--.-3!  hollow  braided  conductor  of  the  type  used 
for  shielded  cable  has  essentially  the  same  high  frequency  resis- 
tance as  a solid  conductor,  but  with  considerably  greater  flexibility. 

1 

Ribbon  Conductors--A  very  flexible  conductor  can  be  made  by 
laminating  a copper  ribbon  between  thin  sheets  of  plastic.  The 
desired  electrical  and  mechanical  properties  might  be  obtained 
with  a 1/8  inch  wide,  0.001  inch  thick  copper  ribbon  sandwiched 
between  1/4  inch  wide  polyethylene  strips.  This  type  of  conduc- 
tor would  tend  to  lie  nearly  uniformly  near  the  bottom  of  a i inch 
tube.  The  edges  of  the  polyethylene  strip  can  be  deformed 
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periodically  during  manufacture  to  provide  support  at  discrete 
points.  Another  technique  for  providing  periodic  support  points 
is  to  put  an  intentional  twist  in  the  ribbon  conductor  as  the  trans- 
ducer is  fabricated. 

Formed  Conductors--A  helical  inner  conductor  configuration  can 
be  fabricated  by  winding  a relatively  springy  conductor  such  as 
solid  copper  or  copperweld  wire  around  a mandrel  and  then 
stretching  it.  If  the  diameter  of  the  resulting  helix  is  slightly 
smaller  than  the  tubing,  it  will  allow  relative  motion  of  the  con- 
ductors, but  will  limit  the  variations  in  gap  to  the  difference  in 
diameters. 

Insulating  Supports  — Various  forms  of  insulating  support  struc- 
tures, such  as  beads  on  the  inner  conductor  or  deformations  in 
the  insulation,  should  also  be  considered. 

Eccentric  Hollow  Dielectric--The  nominal  inner-conductor  spacing 
can  be  maintained  by  a dieletric  with  an  eccentric  hole,  as  shown 
in  Figure  26.  If  the  hole  is  made  slightly  larger  than  the  inner 
conductor,  the  conductor  is  free  to  move  but  its  maximum  devi- 
ation is  restricted  to  an  acceptable  value.  An  inexpensive,  low- 
loss  dielectric  material  such  as  solid  or  cellular  polyethylene  will 
not  significantly  increase  the  attenuation  as  compared  to  an  air- 
insulated  transducer.  Any  of  the  conductor  configurations  described 
above  can  be  used.  This  geometry  is  attractive  because  it  closely 
resembles  the  basic  wire-in-tube  configuration  and  should  be 
capable  of  low-cost  fabrication. 
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Figure  26.  Eccentric  Hollow  Dielectric  Transduc 
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It  is  difficult  to  predict  analytically  how  these  conductor  configurations  will 
respond  to  a mechanical  disturbance.  We  recommend  that  short  segments 
(10  meters  or  less)  of  several  transducer  designs  be  fabricated  for  experi- 
mental determination  of  their  mechanical  and  electrical  characteristics. 

.A  larger-scale  feasibility  demonstration  should  then  be  conducted  with  a 
selected  transducer  design. 
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APPENDIX  A 
SOIL  DISPLACEMENT  COMPUTATION 

In  estimating  the  soil  deformations  caused  by  an  intruder,  we  made  use  of  a 
computer  model  which  had  been  developed  at  Honeywell  prior  to  this  study. 
The  approach  taken  in  computing  soil  displacements  is  to  assume  a quasi- 
static model.  This  model  assumes,  in  effect,  that  input  load  frequencies  are 
relatively  small  so  that  soil  displacements  are  proportional  to  and  in  phase 
with  the  applied  loads.  This  assumption  is  valid  where  the  wavelength  for 
displacement  propagation  is  much  larger  than  that  of  the  ranges  for  detection 
under  consideration.  Most  soils  exhibit  wave  velocities  well  in  excess  of  100 
feet  per  second.  Therefore,  the  one-hertz  input  frequency  results  in  a 100- 
foot  wavelength.  Since  ranges  for  detection  are  for  human  intruders,  well 
under  10  feet,  the  quasi-static  assumption  is  valid  for  that  type  of  traffic. 

As  consideration  is  given  to  higher  frequencies  and  larger  ranges  for  detec- 
tion, i.e.,  as  with  vehicular  traffic,  the  quasi-static  assumption  becomes 
invalid.  One  must  then  resort  to  analysis  of  wave  propagation  within  the  soil. 

ANALYTICAL  SOLUTION  FOR  STRAIN 

The  solution  will  be  worked  out  for  the  case  of  cylindrical  coordinates 
assuming  axisymmetry.  The  problem  is  thus  cast  in  two  dimensions,  r 
and  z,  as  labeled  in  Figure  A-l.  The  load  produced  by  the  intruder  must 
be  assumed  to  act  vertically  at  r = 0. 

For  the  case  of  uniform  elastic  soil  properties  and  a concentrated  vertical 
load  P at  r = 0,  a closed-form  solution  may  be  found  in  several  texts.  The 
solutions  for  radial  and  vertical  displacements  are: 
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where 

E = Young' s modulus;  v = Poisson' s ratio. 

The  assumptions  of  uniform  soil  properties  and  a point  load  limit  the  applic- 
ability of  the  above  equations.  The  point-load  assumption  is.  according  to 
Saint- Venants'  principle,  valid  for  distributed  loads  at  distances  from  the 
load  large  as  compared  with  the  dimensions  of  area  of  application.  Thus, 
for  soil  displacements  in  the  immediate  vicinity  of  a footprint.  Equations 
(A-l)  and  (A-2)  would  give  imprecise  results.  The  assumption  of  uniform  soil 
properties  does  not  conform  to  usual  conditions  of  terrains.  Soil  properties 
vary  considerably  with  depth  either  because  of  harder  strata  existing  at 
lower  depth  or  because  of  compaction  caused  by  overpressure. 


NUMERICAL  SOLUTION  FOR  DISPLACEMENT 

For  the  case  of  variable  soil  properties,  and  for  displacements  in  the  vicinity 
of  a distributed  load,  one  must  resort  to  a numerical  solution.  Such  a solu- 
tion has  been  derived  and  programmed  in  Fortran  IV  on  the  Honeywell  6080 
computer. 

The  basic  methodology  is  to  write  numerical  approximations  for  the  equations 
of  elasticity.  We  are,  of  course,  assuming  the  soil  behaves  according  to 
Hook's  law  (i.  e.  , stress  is  proportional  to  strain).  For  certain  soils. 
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sufficient  water  content  may  produce  a more  plastic -type  behavior  which 
cannot  be  well  represented  by  this  analysis. 


For  cylindrical  coordinates  and  axisymmetry  as  illustrated  in  Figure  A-l, 
the  equations  of  elasticity  with  variable  soil  properties  are: 


ri 

I du.  u | , d2u  1 

l 2 

Idr  * r r rfJ 

dr 


dz2 


a 


+ G)  [ 


d u 
dr  dz 


+;igl|u  + |w 


= 0 


, dz  V dz  dr 

dr  dz 

,[ld\|fdu  u|  , { 0 dG  l d'w  ' = 0 

\\  dzlldr  + r J Mdz  + 2dz  Jd7  ' 


(A-3) 


(A-4) 


where: 


\ 


2E 

(1  + 2)  ( 1 - 2v) 


G 


E 

2(1  + v) 


The  terms  highlighted  by  dashed  lines  account  for  property  variations  in 
the  vertical  direction. 

The  above  equations  have  been  expressed  in  finite  difference  form  and 
programmed  for  solution  using  a relaxation  technique. 

The  program  is  set  up  with  a graded  grid  system:  that  is,  as  the  distance 
from  the  point  of  application  increases,  the  size  of  the  grid  also  increases 
as  shown  in  Figure  A-2.  The  basic  program  is  set  up  for  five  gradations 
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so  that  displacements  are  computed  to  a range  of  128  inches.  A 200-pound 
load  is  assumed  to  act  uniformly  over  a circular  disk  with  a radius  of  four 
inches.  Printouts  of  displacements  are  made  for  each  16  by  16  grid  so  that 
a high-resolution  view  of  displacements  can  be  obtained  in  the  vicinity  of 
the  load.  Figures  A -3  and  A -4  are  computed  values  of  w and  u respectively 
out  to  a radius  of  eight  inches;  Figures  A-5  and  A-6  are  for  a grid  16  times 


0.5  o.  -47.  1 -73.8  -93.2  -109.9  -124.3  -136.3  -145.9  -150.9  -153.2  -152.5  -148.5  -143.7  -138. 


o 


o 


Figure  A-4.  Horizontal  Soil  Displacement,  u (nin.  ),  0.5  Inch  Resolution, 
G = 5000  psi,  v=  0.35 
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Figure  A-6.  Horizontal  Soil  Displacement,  u (u  in.  ),  8-Inch  Resolution, 
G = 5000  psi,  v = 0.  35 


APPENDIX  B 


j 


CHARACTERISTIC  IMPEDANCE  OF  AN  ECCENTRIC  COAXIAL  CABLE 
WITH  A NONHOMOGENEOUS  DIELECTRIC 


THE  INVERSION  TRANSFORMATION  OF  A CIRCLE 

Consider  the  properties  of  a regular  function  of  the  type 

Z = f(t)  = x(u»  v)  + jy(ii,  v)  (B-l) 

which  defines  a complex  variable  Z = x + jy  as  some  function  of  another 
complex  variable  t = u + jv.  A particular  value,  t\  can  be  represented  by 
a point  in  the  Argand  diagram  of  t.  Through  equation  (B-l)  some  particu- 
lar value  Z'  corresponds  to  t'  and  it  may  be  represented  by  a point  in  the 
Argand  diagram  of  Z.  Further,  there  is  a similar  correspondence  for  a 
succession  of  pairs  of  points  in  the  t and  Z planes  so  that  to  some  curve  t't" 
there  corresponds  a curve  Z'Z"  which  is  said  to  have  been  transformed 
from  t't"  by  equation  (B-l). 

Consider  the  inversion  transformation 

Z = 1/t  (B-2) 

The  equation  of  a circle,  center  (a.o),  and  radius  r in  the  t plane, 

r = ! t-a  | (B-3) 


Substituting  (B-2)  into  (B-3), 


|_J— 

(x+jy) 


,(l-ax)  - ja} 
x+jy 


(B-4) 


Therefore, 


2 2 2 

2 (1-ax)  + a y 


2 2 
x + y 


(B-5) 


which  may  be  rearranged  as 


x’“2 2 +y 


2 2 
a - r 


(B-6) 


2 2 2 2 

Equation  (B-6)  represents  a circle,  radius  r/(a"-r‘"),  center  ((a/(a“-r  )),  0), 
in  the  Z plane,  mapped  by  equation  (B-2)  from  the  circle,  radius  r,  center 
(a,  o),  in  the  t plane. 


THE  IMPEDANCE  OF  AN  ECCENTRIC  COAXIAL  CABLE 
WITH  A NONHOMOGENEOUS  DIELECTRIC  FILLING 


Consider  the  eccentric  coaxial  cable  shown  in  Figure  B-l.  From  equation 
(B-6),  it  is  seen  that  the  smallest  circle,  center  (d,  0),  and  radius  R^  in  the 
t plane  transforms  into  a circle  in  the  Z plane  with  center  coordinates 
((d/(d2-R32)),  0)  and  radius 


(B-7) 
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DIELECTRIC  WITH  DIELECTRIC 
CONSTANT  e. 


CONDUCTOR 


DIELECTRIC  WITH  DIELECTRIC 
CONSTANT  c„ 


Figure  B-l.  Eccentric  Coaxial  Cable  with 
Nonhomogeneous  Dielectric 


Similarly,  the  circle  with  center  (d,  0)  and  radius  R in  the  t plane  trans 

„ ^ 9 2 

form  into  a circle  in  the  Z plane  with  center  coordinates  ((d/(d“-R9  )),  0) 


r = 

2 2 9 
- d2-p- 


.Also,  the  circle  with  center  ((d-D),0)  and  radius  R.  in  the  t plane  trans- 

2 

forms  into  a circle  in  the  Z plane  w'ith  center  coordinates  ((d-D)/(d-D)  - 
2 

R.“  , 0)  and  radius 


For  the  conductors  to  be  concentric  in  the  Z plane,  it  is  necessary  that 


d _ d-D 
d2  - R2  (d-D)2  - R2 


(B-1C 


2D 


Equation  B-10  fixes  the  necessary  position  of  the  origin  of  inversion  in  the 
Z plane  by  defining  d as 


D2  + R2  - Fj2  2 


\ 


2D 


I 


R. 


(B-ll 


The  value  of  d making  the  r.(i  = 1, 2,  3)  positive  must  be  taken 


y 


DIELECTRIC  WITH 
DIELECTRIC  CONSTANT  e2 


Figure  B-2.  Transformed  Cable 


The  configuration  of  Figure  B-2  can  be  broken  up  into  infinitesimal  wedges 
as  shown  in  Figure  B-3. 


TYPICAL  WEDGE 
(EXAGGERATED) 


m(e) 

V 


d-R^  d-R^ 


Figure  B-3.  Wedge  Configuration 


The  approximate  capacitance  per  unit  length  of  the  cable  is 


2£°  ) \ 1/£1  Xn  ^ . 1/C2  ^ 

0 \ m (0)  , 


(B-12) 


The  inductance  per  unit  length  of  the  cable  is 


^0  ri 
l*  — In  — 

2TT  r„ 


(B-13) 


The  characteristic  impedance  of  the  cable  is 

Z - VT"  (B-14) 

c" 

where  C and  l are  respectively  defined  in  equations  (B-12)  and  (B-13). 

Using  this  technique,  we  computed  the  characteristic  impedance  of  a trans- 
ducer consisting  of  a 0.018-inch  diameter  wire  in  a 0.  19-inch  ID  tube.  The 
outer  diameter  of  the  insulation  was  0.03  inches  and  its  relative  permea- 
bility was  2.  18.  With  the  insulated  wire  touching  the  outer  conductor,  the 
characteristic  impedance  is  47  ohms.  The  same  conductor  spacing  would 
give  a characteristic  impedance  of  62  ohms  for  an  air-insulated  transducer, 

and  42  ohms  for  a uniform  dielectric  with  e =2.18. 

r 

When  there  is  a finite  air  gap,  the  effect  of  a thin  dielectric  layer  on  the 
characteristic  impedance  is  considerably  reduced.  The  air-dielectric 
approximation  for  the  transducer  impedance  then  becomes  reasonably 
accurate.  * 
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APPENDIX  C 

REFLECTIONS  FROM  ABRUPT  AND  GRADUAL  DISTURBANCES 
IN  A LONG  TRANSMISSION  LINE 


This  appendix  describes  the  equations  used  in  predicting  reflected  signal 
waveshapes  for  a uniform,  lossy  transducer.  Our  computations  were  simpli- 
fied by  using  a double-step  approximation  for  the  line  displacement.  We 
include  in  this  appendix  an  independent  computation  of  the  signal  waveshape 
for  a case  in  which  the  impedance  variation  is  described  as  an  exponential 
function.  Results  of  the  two  computation  methods  are  nearly  identical  for 
long  transducers,  where  frequency-dependent  losses  dictate  the  signal 
waveforms. 


STEP-FUNCTION  RESPONSE  OF  A TRANSMISSION  LINE 


According  to  the  theory  developed  by  Wigington  and  Nahman  (Reference  1). 
the  response  of  a lossy  transmission  line  to  a step-function  input  is 


v(t)  = CERF  y 


(C-l) 


where  t is  the  time  after  arrival  of  the  step  and  T is  defined  by 


T = 


2 2 
or  L 


4 Z 


(C-2) 


L is  the  length  of  the  transmission  line  and  Z is  its  characteristic  imped- 

o r 

ance.  The  attenuation  per  unit  length  a is 
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a = -LI  I—  + — 

2n\-"l  K a2 


and  a9  are  the  conductor  radii,  and  for  copper  conductors 


(i  = p.  = 4TT  x 10  Henries/meter 
o 

7 

a = 5.  8 x 10  mhos/meter 


CERF(x),  the  complementary  error  function  of  x,  can  be  computed  very 
accurately  from 


CERF(x)  = 


, 2 3 41  4 

1 + a^x  + a9x  + a3x  + a^x 


(C-4) 


a = 0.  278393 
a9  = 0.  230389 
a = 0.000972 

i) 

a . = 0.  078108 
4 


REFLECTIONS  FROM  A DOUBLE-STEP  DISTURBANCE 
IN  LINE  IMPEDANCE  (METHOD  1) 


When  an  intruder  crosses  the  transducer,  the  impedance  of  a short  section 
is  affected.  In  our  analysis,  we  assumed  that  the  disturbance  can  be  approx- 
imated by  step-function  changes  in  the  line  impedance  (Figure  C-l).  Let  us 
denote  the  line  impedance  of  the  disturbed  section  as  Z ',  and  the  length 
Ax.  We  now  want  to  calculate  the  reflected  wave.  This  is,  in  fact,  a 
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multiple  reflection  problem.  However,  we  can  make  a "one  pass"  approx- 
imation by  realizing  the  fact  that  Zq  and  Z ' are  almost  equal;  therefore, 
the  transmission  is  nearly  complete.  Let  us  assume  that  an  incident  wave 
v(t)  arrives  at  junction  1.  The  reflection  coefficient  is: 

Z ' + Z 
- o o 

P1  Z + Z ' (C-5) 

o o 

The  transmission  coefficient  is: 

2Z  ' 
o 

T1  " Z + Z ' 
o o 


(C-6) 


This  transmitted  wave  propagates  forward  and  is  reflected  at  junction  2. 
The  reflection  coefficient  is: 


Z ' - Z 
o o 

P2  = Z + Z ' 


(C-7) 


This  reflected  wave  propagates  backward  and  is  transmitted  at  junction  1, 


where  the  transmission  coefficient  is: 


T i — 

1 Z + Z ' 
o o 


Now  this  wave  has  travelled  a distance  of  2Ax,  and  the  time  spent  is 


(C-8) 


(C-9) 


Combining  Equations  (C-6),  (C-7)  and  (C-8),  the  reflected  wave  has  a mag- 
nitude 


v (t  + At)  = T1T2T2,v(t) 


(C-10) 


However,  at  t + At,  the  r.ewly  arrived  wave  v(t  + At)  is  also  reflected,  the 
reflected  wave  magnitude  is: 


v"(t  + At)  = p jv(t  + At) 


(C-ll) 


Combining  Equations  (C-10)  and  (C-ll),  and  using  p 2 ~ ~P^»  total 
reflected  wave  magnitude  is 


+ At)  = p^vft  + At)  - r1T2'v(t)] 
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Equation  (C-12)  becomes 


vr(t  + At)  « pjMt  + At)  - v(t)]  (C-15) 

This  is  equivalent  to 

vr(t)  * pj[v(t)  - v(t  - At)]  (C-16) 

We  can  consider  this  as  a superposition  of  two  separate  reflected  signals 
displaced  in  time.  Each  signal  is  a linearly-scaled  version  of  the  incident 
wave  v(t).  As  these  signals  propagate  back  toward  the  generator,  they 
continue  to  undergo  the  same  degradation  in  risetime  as  the  forw'ard-going 
incident  wave.  We  can  compute  the  actual  waveshape  of  the  superposed 
reflected  signals  by  using  the  expression  previously  developed  for  v(t) 
evaluated  at  a distance  L equal  to  twice  the  distance  from  the  generator  to 
the  disturbance. 

REFLECTIONS  FROM  AN  EXPONENTIAL  VARIATION 
IN  LINE  IMPEDANCE  (METHOD  2) 

As  a check  on  the  validity  of  the  double-step  approximation  for  the  trans- 
ducer impedance  variation,  we  computed  the  signal  waveshape  from  an 
exponential  disturbance.  The  expression  for  the  variation  in  impedance 
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was  taken  from  a cable  sensor  analysis  performed  by  Frazier  and  Shiau 
(Reference  2).  They  expressed  the  impedance  of  the  disturbed  section  as 

Z (x)  = Z exp  [-k(x  - L)/Ax],  L < x < L + ^ (C-17) 

O O u 

Z (x)  = Z exp  {-k[l  - (x  - L)/Ax]},  L + <x<L  + Ax  (C- 18) 
o o ^ 

The  parameter  k is  determined  from  the  impedance  at  the  center  of  the 

disturbed  section  Z ' by 
o J 

Z 

k = 2*n  ^-2.  (C-19) 

o 

Note:  For  small  variations  in  impedance,  the  expression  for  Zq(x)  is 
nearly  a triangular  function. 

For  time  domain  reflectometry  purposes,  the  reflected  waveform,  v^lt),  in 
a uniform  transducer  with  an  exponential  disturbance  when  interrogated  by 
a unit  step  voltage  is  given  by  the  expression 

2 

vr(t)  = ^n^r  [u(t)  * F_1  {[exp  (4vL-vAx)-exp(4vL)]2/v}]  (C-20) 

o 


where 


L = the  distance  between  the  reflectometer  and  the  disturbance 

Ax  = the  length  of  the  disturbance 

ZQ  = the  characteristic  impedance  of  the  transducer 

ZQ'  = minimum  impedance  along  the  exponential  disturbance 

Y = propagation  coefficient  of  the  transducer 

u(t)  denotes  a unit  step  function 

F * is  the  inverse  Fourier  transform 
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The  propagation  coefficient,  y»  of  the  transducer  is  given  by  the  expression 


V (jw)  = ^(R  + jwi  )/jwc 


where  l - inductance  per  unit  length  of  the  transducer 
c = capacitance  per  unit  length  of  the  transducer 
R = k(jw)* 

k » JL 1/2  l±+±\ 


r^  = inner  conductor  radius 
r 2 = outer  conductor  radius 
a = conductivity  of  the  inner  conductor 
p,Q=  permeability  of  free  space 


COMPARISON  OF  RESULTS 

Figure  C-2  shows  a reflected  signal  waveform  computed  by  the  first  method, 
assuming  a step  discontinuity  in  the  line  impedance.  The  disturbance  is 
located  at  one  km  from  the  generator  end  of  a transducer  with  a 1 = 4.  2 mm 
and  a^  = 15  mm.  The  gap  between  conductors  is  0.06  mm,  except  for  the 
1 -meter-long  disturbed  section  where  the  gap  is  0.  05  mm.  A 1-volt  step 

is  used  as  the  excitation  signal.  For  this  case,  the  peak  ampli  ude  of  the 

-4 

reflected  pulse  is  -4.  1x10  volts,  and  the  half-amplitude  width  is  about 
-7 

3 x 10  seconds. 


For  comparison,  we  used  the  second  method  to  compute  the  reflected  signal 
from  an  exponential  impedance  variation  on  an  identical  transducer.  The 
total  disturbance  length  was  again  chosen  to  be  one  meter.  We  set  the 
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impedance  Z ' in  the  center  of  the  disturbed  section  to  be  equal  to  Z ' 

° o 

for  the  previous  case. 

A 512-point  Fast  Fourier  transform  and  discrete  convolution  with  the  unit 
step  input  were  performed  on  a SDS  9300  computer.  Results  of  the  compu- 
tation over  the  time  interval  of  interest  are  shown  in  Figure  C-3.  The 

-4 

peak  amplitude  of  -2.2  x 10  volts  is  about  half  the  amplitude  for  the 
double-step  case,  but  the  pulse  width  is  essentially  unchanged. 

We  would  obtain  almost  identical  results  with  Method  1 if  we  replaced  the 

impedance  Z ' of  the  one-meter  disturbed  section  bv  (Z  ' + Z )/2.  Inte- 
o o o 

grating  the  exponential  expression  for  the  variation  in  transducer  imped- 
ance over  the  disturbance  length,  we  find  that  the  average  value  is  in  fact 

very  close  to  (Z  ' + Z )/2. 

o o 

This  comparison  of  numerical  results  does  not  constitute  a proof  of  the 
validity  of  Method  1.  However,  we  feel  that  the  method  provides  accurate 
results  for  long  transducers,  where  the  response  in  the  time  domain  is 
limited  by  the  frequency  response  of  the  transducer.  Under  these  condi- 
tions, the  amplitude  of  the  received  signal  from  a disturbance  will  depend 
upon  the  integral  of  the  impedance  variation  over  the  disturbed  length. 


I 
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Figure  C-3.  Signal  Waveform  Computed  by  Method  2 
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METRIC  SYSTEM 


BASK  UNITS: 

Quantity 


Unit 


length 

mui 

time 

electric  currant 
thermodynamic  temperature 
amount  of  subetanca 
luminous  intensity 

SUPPLEMENTARY  UNITS: 

plane  angle 
solid  angle 

DERIVED  UNITS: 
Acceleration 

activity  (of  a radioactive  source) 
angular  acceleration 
angular  velocity 

area 

density 

electric  capacitance 

electrical  conductance 

electric  field  strength 

electric  inductance 

electnc  potential  difference 

electric  resistance 

electromotive  force 

energy 

entropy 

force 

frequency 

illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

quantity  of  electricity 
quantity  of  heat 
radiant  intensity 
specific  heat 
stress 

thermal  conductivity 
velocity 

viscosity,  dynamic 

viscosity,  kinematic 

voltage 

volume 

wavenumber 

work 


metre 

kilogram 

second 

ampere 

kelvin 

mole 

candela 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  second  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

loule 

joule  per  kelvin 

newlon 

hertz 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

coulomb 

loule 

watt  per  steradian 

loule  per  kilogram-kelvin 

pascal 

watt  per  metre-kelvin 
metre  per  second 
pascal-second 
square  metre  per  second 
volt 

cubic  metre 
reciprocal  metre 
loule 


SI  Symbol  Formula 

m 

kg 

a 

A 

K 

mol 

cd 


rad 

sr 


m/s 

(disintegretionVs 

rsd/s 

rad/s 

m 

kgim 

F 

A-s/V 

S 

A/V 

Vm 

H 

V-s/A 

V 

W/A 

ViA 

V 

WiA 

| 

N-m 

1 K 

N 

kg- ova 

Hz 

(cycle  Vs 

lx 

im/m 

cdim 

Im 

cd-sr 

Am 

Wb 

V-s 

T 

Wbrtn 

A 

W 

1/* 

Pa 

Nm 

C 

1 

A-s 

N-m 

War 

I'kg-K 

Pa 

N'm 

Wm-K 

m/s 

Pa-s 

m/s 

V 

W/A 

m 

i 

(waveym 

N-m 

SI  PREFIXES: 


Multiplication  Factors 


Prefix  SI  Symbol 


1 000  000  000  000  - 

in11 

tera 

T 

G 

M 

1 000  000  000  - 

10* 

*'*• 

t 000  000  - 

10* 

mega 

t 000  * 

10* 

kilo 

k 

too  • 

101 

hecto* 

h 

to  • 

10' 

deka* 

da 

0 1 - 

10-' 

decl* 

d 

0 01  • 

10-* 

centi* 

c 

0 001  » 

10*‘ 

milll 

m 

0 000  001  • 

I0-* 

micro 

0 000  IKK)  001  • 

10-* 

nano 

n 

0 000  000  000  001  • 

10-'1 

pico 

p 

0 000  000  000  00(1  001  « 

10-'* 

femto 

V 

0 000  000  000  000  IKK)  001  « 
To  be  avoided  where  possible. 

lO-'F 

alto 

• 
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